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Abstract 


It is the purpose of this research to examine the engineering process required to 
adequately predict the total resistance of unique hullforms. A companion of the significant 
physical effects present in monohull and advanced hullforms testing is made. The histoncallv 
accepted method of scaling using Froude’s Hypothesis is studied and its underlying 
assumptions are considered in detail. An overall design philosophy is created for achieving 
resistance minimization through a range of operating speeds. This discussion is applied to a 
specific hullform called SLICE which is an adaptation of the Small Waterplane Area Twin 
Hull (SWATH) concept. 


To further develop the body of available SLICE resistance data. three models were 
built and tested at the US Naval Academy's 380 feet towing tank in Annapolis Maryland. The 
results of these tests 1s reported on and a video demonstration of their performance 1s 
provided as part of this report. 


Finally, a numencal method approach to the solution of Mitchell's integral is presented 
as a possible means of calculating SLICE wave making resistance. This is considered as an 
intermediate step tool in design optimization, useful after preliminary model tests are 
completed and before final designs are decided for powered model testing. 


Thesis Supervisor: Dr. Alan J. Brown 
Title: Professor of Naval Architecture and Manne Engineering 
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Chapter | Introduction 


The majority of the world's commerce is moved by monohull displacement 
ships. Ironically, in this age of supersonic air travel, guaranteed overnight mail 
delivery and fibre optic data transfer at the speed of light, modern displacement ships 
are capable of speeds no faster than those built fifty years ago. The apparent lack of 
progress in high speed ocean travel stems from the phvsical realities of the interaction 
between the displacement monohull form and the water in which tt travels. 

The total resistance curve for a ship 1s generally plotted as the Total Drag 


Coefficient. C,; versus Froude Number. . These are expressed as: 


F=V/STQLY cy 
Ry=(5 PVC, (1-2) 


where L is some characteristic length. usually length of the waterline of the body, p 1s 
the density of the water. V is the ship's velocity and S 1s the Wetted Surface Area. 
Total resistance 1s generally recognized as the additive effects of friction. wave 
making and spray generation of both the main bodv and its appendages. The extent to 
which each of these components effects the total resistance 1s. however. the subject of 
ongoing debate even though it has been over one hundred years since Froude gave us 
this theory. For a monohull displacement ship, the wave making resistance portion of 
the total resistance curve has the characteristic shape shown in Figure (1-1). Since 
resistance is obtained by multiplying the coefficient shown in Figure (1-1) by V’, the 
total resistance curve becomes nearly vertical and practically insurmountable by 


conventional hull forms at the speeds shown. 






RESISTANCE COEFFICIENT 


SLUICE (20 kts) 40 ® 


0 Q4 0.5 19 13 
FROUDE No. 2 Vv VoL 


ere ee ee FE=s 


Figure I-1: Characteristic Wave Making Resistance Curve 
(from LNISC presentation) 


For Naval Architects. the search is ongoing in advanced hull form design for a 
practical method to “get over the hump”. fost advanced hull forms do not resemble 
the classic displacement monohull shape on which the vast majontv of Naval 
Architecture experience is based. The purpose of this research is to examine the 
underlying assumptions regarding commonly accepted full scale resistance prediction 
methods in an effort to understand the important features of design optimization and 
draw some conclusions regarding resistance minimization methods . The Lockheed 
Missile and Space Company in conjunction with Pacific Marine (L\{SC,PM) has 
proposed a hull form which is a modification of Small Waterplane Area Twin Hull 
(SWATH) concepts. Their proposal. called SLICE. will be presented to the United 
States Government as a 177 Lton vessel in an Advanced Technology Demonstration 
(ATD) in June 1996. A sketch of this hullform is shown in Figure 1-2. The SLICE 
vessel will be the subject hullform for the research conducted as part of this paper. 


SLICE is protected bv Lockheed under patent law. 
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Figure 1-2: Artist's Concept of SLICE from an unpublished LMSC, 
Pacific Marine. ONR Article. 


1.1 The SLICE Concept 


High speed is not an inherent quality of SWATH ships. Although there may 
be some reduction of the wave making resistance due to smaller waterplane area. there 
is also an increase in the total wetted surtace which increases frictional resistance. In 
SWATH ships. high speed is a design teature. When compared to an equivalent 
displacement monohull. SWATH design must be considered as a trade-off between a 
definite increase in frictional resistance and the potential for a decrease in wave 
making resistance. The design question becomes how to ensure that the decrease in 
Wave making resistance is greater than the gain in frictional resistance. 

LMSC's SLICE concept differs trom traditional SWATH design in that the 
demi-hull is replaced bv two shorter pods. Each pod is attached to the main box 
structure by a single strut. It is through this design feature that LMSC intends to 
minimize the impact of wave making resistance. Bv designing the vessel with an 


unusually short characteristic length. the hump of the wave making curve will in 


~ 
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theory be driven to a very low speed regime. Since resistance is a function of V*, any 
left shift of the curve could have a substantial minimizing effect on the total 
resistance associated with the hump. It is further theonzed that in the high speed 
regime. total resistance will be dominated by fnctional effects and therefore design 
efforts focused on minimizing the wetted surface of the hullform will vield the best 
result. LMSC has chosen the length of the submerged pod as the characteristic length 
critical to design considerations. It should be emphasized that this hullform does not 
eliminate the hump effect characteristic of the wave making resistance curve. It does 
not specifically cause wave cancellation effects nor does it attempt to reduce viscous 
drag effects. Rather, the apparent thrust of the LMSC theory 1s to stagger the peak 
effects of wave and skin resistance. thereby reducing the overall resistance value at 


specific design speeds. 


1.2 LMSC Model Test Synopsis 

LMSC has conducted two different sets of model tests at the US Navy's 
Surface Warfare Center, Carderock Division. Bethesda, Marvland (NSWC). In the 
first series of tests. called the HM&E tests, several different configurations of struts, 
pods and foils were used to make a determination of the overall best configuration to 
use as the baseline for the production vessel. These tests were conducted in the Spring 
1994. 

The second set of tests. also conducted at NSWC. were used to obtain detailed 
resistance and control surface performance data for the baseline production vessel. 
These tests. conducted in August 1994 and continued in December 1994 were 
observed by several MIT students. The model that was tested was a 1/8th scale of the 
production vessel. It's construction was fiberglass bodies mounted to an aluminum 
box frame. For all the tests conducted, the model was fixed in pitch and heave and 
data was collected for forces and moments about the three principle axes. There were 
two major groupings of data collection runs performed. The first group was called the 
dynamic control tests and was primarily used to judge the required effectiveness of 


control surfaces to hold the ship at a constant draft with an even tnm. Each mun had a 
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slightly different model configuration with canards, stabilizers and rudders tested in 
different locations and with varving angles of attack. Each configuration was tested a 
single time unless there was some specific problem noted in the conduct of the test. 
The second major grouping of tests consisted of rather major modifications to the 
model configuration in an attempt to improve the design. Haunches on top of the stnuts 
were removed to allow for deeper strut length, wedges located at the top of the struts 
were added and removed to experiment with better surface piercing shapes, the overall 
draft of the model was adjusted to force a deeper submergence of the rather bluff pods 
and finally. a new model was tested in which the longitudinal position of the struts 
relative to the pods was altered to modify the flow over the pods. 

Within these two major test categories. the test procedure remained consistent. 
The model remained fixed in pitch and heave for all tests. Two different speed 
regimes were tested. a lower speed run consisted of twelve speed increments ranging 
from 3.2- 8.8 feet per second (fps) corresponding to full scale speeds of 9-25 knots. A 
high speed run used seven increments from 5.7-12.4 fps which corresponds to full 
scale speeds of 10-35 knots. All speed increments were fit into a single pass along a 
tow tank length of approximately 1200 feet. All of the dynamic control tests were 
conducted only in the lower speed regime. the major model modification tests were 


conducted in both speed regimes. 


1.3 Research Motivation 

The US Navy has potential applications for relatively small (200-700 Lton) 
high endurance craft that embody the seakeeping attributes of SWATH. The Navy's 
interest in the SLICE program is demonstrated by the funding the Navy has committed 
to the ATD that LMSC will perform in June 1996. It was LMSC/PM's entrepreneurial 
spirit that initially brought the SLICE concept to the attention of the Navy but this 
same entrepreneurial spirit has the potential to limit the usefulness of the SLICE 
design as well. LMSC.PM's allocated timetable from concept development to final 
product delivery is just over two years, very short by Navy shipbuilding standards. As 


a result LMSC must focus on a relatively narrow set of engineering problems that 


specifically address the needs of the 177 Lton ATD vessel. Broader based questions, 
such at an in-depth understanding of the effects of body interactions. how these 
interactions vary with modifications to the basic dimensions, the practicality of scaling 
the SLICE concept into a different size regime (1500 Lton. 5000 Lton. etc.) and a 
relative comparison of the merits of a SLICE ship as compared to a SWATH or 
Monohull design for a given mission are all facets of the development of this concept 
that are not specifically addressed by LMSC in the ATD. 

This thesis will attempt to add to the developing body of SLICE knowledge bv 
focusing on the resistance and powering aspects of the SLICE vessel in an 
academically mgorous manner. In so doing, three benefits will be realized: 

1) Independent research into this hull form will give both the US Navy and 
LMSC a second hydrodynamic viewpoint of the 177 Lton vessel prior to the ATD. 

2) Development of a sufficient body of data will aid optumization of the 177 
Lton hull form. 

3) A broader understanding of the hydrodynamic features of the hull form 
will facilitate deciding the merits of pursuing designs in other displacements more 


practical to the US Navy. 


1.4 Approach 
To achieve the objectives outlined above. there are several questions raised by 
the LMSC efforts that require answers. These questions and the method that will be 


pursued for answering them is outlined in this section. 


ire ORY 

LMSC theorizes that the appropriate characteristic length for SLICE is the 
length of a single pod. In conventional designs. is based on the length of the 
waterline (LWL). Representative characteristic lengths for SLICE would seem to be 
bounded on the low end by the length of a single strut and on the high end by the 
length between the leading edge of the forward strut to the trailing edge of the aft 


strut. LMSC's choice of an underwater body for the charactenstic length is not an 
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obvious one. Understanding the appropriate charactenstic length for the SLICE 
configuration 1s necessary to understand the magnitude of savings that can be expected 
over competing hull forms. 

The interaction between bodies is the second area for further study. There are 
two specific situations that require detailed analysis. The first is the effect a strut and 
its attached pod have on one another. the second is the effect that a strut/pod 
combination have on any of the other three strut/pod combinations that form the 
underwater body. Understanding these effects is necessary for design optimization and 
also deciding the practical upper limit of displacement for a SLICE ship. 

Finally, the LMSC proposed design places the forward struts inboard of the aft 
struts as illustrated in Figure 2. This arrangement may cause some significant 


outward lateral forces on the struts that must be considered in the structural design. 


TESTING 

There were several procedural questions raised by the model tests conducted by 
LMSC in August. Keeping the model fixed in pitch and heave was the most 
significant issue. Although it was intentionally done to facilitate the dynamic control 
analvsis, it also left open significant questions regarding how the hull form behaves at 
different speeds. Considering the recognized sensitivity that SWATH type hulls have 
in tons per inch immersion (TPI) and the moment to tnm one inch (MT1) this may be 
significant. 

The short duration of each test run speed increment, the lack of repetition of 
tests. and the truncation of most data collection runs at the full scale equivalent of 25 
knots for a ship required to achieve 30 knots raises questions regarding the 
thoroughness of the results. Although this is understandable considering the expense of 
model testing and the limited time available, more testing must be completed before 


the SLICE body of data can be utilized to make design decisions. 


OPTIMIZING 


There are several striking features incorporated in the LMSC design. The most 


stnking is the bluffness of the leading sections of the pods and struts, especially for a 
high speed craft. It would appear that minimizing the wetted surface to displacement 
ratio was LMSC's dnving concern in hullform selection. Although this approach 
minimizes frictional resistance, it also introduces significant form drag. With a pod 
length to diameter ratio (L/D) of 4.2, SLICE cannot be thought of as a slender body. 
This introduces another concern in evaluating SLICE performance. The vast majority 
of traditional resistance analysis has slender body theory incorporated to some extent. 
Any evaluation process undertaken must be carefully scrutinized to ensure that the 
SLICE hullform does not violate this imbedded assumption. A design approach which 
more thoroughly addresses the entire speed spectrum 1s preferable even with a specific 
mission speed identified. All speeds from zero through the mission speed must be 
attainable so some design consideration must be given to the entire speed spectrum. 
With eight component bodies ( 4 struts and 4 pods), wave cancellation effects through 


body interaction should be a design focus. 


RESEARCH METHOD 

In an attempt to answer these questions the following work was accomplished: 

1) A literature review that encompassed one hundred vears of ship resistance 
research was performed. This was accomplished to ensure that the full body of 
assumptions made in this field of engineering 1s understood. 

>) A study of model testing techniques was made to develop an appreciation 
for proper test methods. 

3) Series 58 resistance data was studied for potential pod optimization. 

4) A literature survey was conducted to ascertain the state of the science in 
spray drag predictions. 

5) Three models were built and tested to expand the overall data bank. 

6) A numerical calculation method for wave making resistance based on the 
solution to Mitchell's integral was attempted in an effort to expedite the hullform 


optimization process. 


Chapter 2. Theory 


The difficulty in predicting full scale resistance is in understanding the extent 
to which different physical phenomena contribute to total resistance. Although effects 
such as boundary layer development. wave making, spray generation and flow 
separation are known to be present as a body moves through the water. methods for 
quantitatively evaluating the contribution of each of these effects to total resistance 
and the extent to which these phenomenon may interact continues to be the subject of 
ongoing debate. 

This chapter studies the two available approaches to resistance prediction, 
model testing and analvsis. The objective is to identify those methods that will be best 
suited for SLICE resistance predictions. Since the vast majoritv of available resistance 
data was developed from displacement monohull analysis. care must be taken to 


ensure that underlying assumptions to these methods are not violated. 


2.1 Model Testing 
The classical approach to model testing 1s well documented. Pnnciples of 

Naval Architecture’ (PNA), and Rawson and Tupper are probably the two most 
recognized authorities outlining this process. The textbook approach 1s definitively 
stated and easv to follow. Nonetheless. the SLICE experience has shown that specific 
model testing methodology and subsequent data reduction can be hotly contested. The 
difficulty lies in deciding if the model tested accurately replicates the physical 
phenomena that will be present in the full scale ship and what to do with the model 
scale data once it 1s obtained. 

As background. the model testing/scaling process outlined by William Froude 


is presented. In this discussion the subscript m refers to the model and the subscript s 


' Comstock. JP. Principles of Naval Architecture, SNAME, New York, 1967. Chapter 7. 
> Rawson.kK J & Tupper.E.C. Basic Ship Theory, 3rd Edition, Vol 2, New York.1983. 
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refers to the full scale ship: 
1) Build a model that is geometrically similar to the full scale ship. This is 
accomplished by choosing a scale factor to apply to all linear dimensions of the ship. 


The scale factor A is then defined as: 


pea (2-1) 


2) Test the model at equivalent speeds. Through dimensional analysis, it 1s 
easily shown that the full scale and model scale speeds are related through equation 


(1-1). Specifically: 


Vel ape (2-2) 


3) Convert the total model resistance measured R,. to a coefficient of total 


tm 


resistance C,. through relationship (1-2). 


tm 
+) Calculate the coefficient of fmetional resistance C,.. The method for 
accomplishing this will be discussed shortly. 
5) A quantity called the coefficient of residuarv resistance C_ 1s obtained bv 


from C,. The reason C,,. is not subjected to the scaling law expressed 


tm tm° tm 


subtracting C 
in (2-2) will be discussed shortly. 
6) The full scale coefficient of residuary resistance C,. is then expressed at the 


same Froude Number as C,. Specifically. 
C =C (2-3) 
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7) The full scale coefficient of frictional resistance C,, is calculated in the 
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same manner as in step 4 and is added to C,, to obtain the full scale coefficient of 
total resistance C... 

8) Full scale total resistance R,. is then found using (1-2). 

This process can be thought of as a system of one equation and one unknown. 
The unknown ts C,,; the equation 1s C,, = C,, + C,,. This method of scaling data is 
useful onlv if the following conditions exist: 

1) A geometrically similar model can be built and tested in compliance 
with the equalities set forth in (2-1) and (2-2). 

2) The model ts tested in a condition that can be considered equivalent 
to the normal operating condition of the full scale ship. 

3) Frictional resistance can be reasonably calculated at both the model 
and full scale size. 

4) All other physical phenomenon present must either obev equality 
(2-3) and therefore be lumped into the residuarv component or be of such little 
concern that ignoring them is inconsequential to the final outcome. 

If the last condition can not be satisfied for a particular phenomenon, then its 
effect must be calculated at both the model and full scale size and algebraically treated 
the same as the friction resistance component. 

Froude's experiments were conducted on slender bodies, generally thought of 
as those with L/D > 10. Residuary resistance was defined as the combined effect of 
wave making resistance and eddy resistance. At low speeds. eddies appear to be a by 
product of the boundary layer and therefore part of frictional resistance. 
Experimentally, Froude demonstrated that eddy resistance did not change 
proportionally with skin friction as parallel midbody was added to a given hull form. 
He concluded that eddy resistance was actually related to the form of the hull. Since 
wave making and eddy resistance are both related to form and are pressure related 


phenomenon. they historically have remained combined as the residuary resistance 


component’. For slender bodies the only significant physical phenomena are friction, 
waves and eddy formation. The process outlined above was developed with slender 


bodies and only these three phenomena in mind. 


2.2 Frictional Resistance 

The problem of calculating the frictional resistance of a ship was first studied 
by considering flat plates. The essential issue focuses on the development of the 
boundary laver along the length of the plate and its transition from a laminar to a 
turbulent regime. Frictional resistance is usually expressed as a coefficient of friction 


in the form of (1-2) and ts generally plotted as a function of the Reynolds Number, §,; 


R= (2-4) 





where L is some characteristic length. usually length of the waterline of the body, v 
1S the kinematic viscosity of the water and V ts the ship's velocitv. The development 
of a laminar flow equation was the work of Blasius (1904). while the turbulent 
solution was formulated independently by Prandtl and von Karman in the 1920's. In 
the range of Reynolds Number between 5x10° and 9x10°, neither the turbulent nor the 
laminar expression vields accurate results. This region of uncertainty 1s known as the 
transition region. Experimentally it has been shown that the transition from laminar to 
turbulent flow occurs gradually across the length of a plate as its speed increases. This 
gives physical meaning to the rather large band of uncertainty between fully laminar 
and fully turbulent flow. The curves of laminar. transition and turbulent plate flow are 
shown in Figure (2-1). 


Model testing through the 1940's attempted to draw parallels between monohull 


: Gilmer,T C and Johnson.B. Introduction to Naval Architecture. Naval Institute Press, Annapolis, 
1982, p.217 


}- 
Cy 


forms and the known flat plate solution. The most accepted flat plate solution to 
frictional resistance adapted to ship form use is the Schoenherr Formula. This flat 
plate formulation is frequently seen in the literature labeled as the ATTC Line. In an 
attempt to account for the effect of hull form on boundary laver formation the 
Intemational Towing Tank Conference (ITTC) adopted an imrerim solution that was 


derived from a correlation of previously obtained model tests. 
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Figure 2-1: Laminar and Turbulent Flow Skin Friction Lines’ 
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This formulation, known as the ITTC 1957 line ts given as: 


Oe 0.075 


a a (2-5) 
* dog) 8-2) 


It should be stressed that this line was developed from a correlation of model data, it 
does not physically represent the behavior of a flat plate and it was not obtained by 
extrapolating flat plate data to ship like forms. Forty years later, the ITTC has not 
replaced this interim solution with a better alternative. 


At the same trme that the [TTC 1957 formulation was published. Hughes 


proposed a method for calculating flat plate frictional resistance yviven bv: 


ce 0.066 
2 eR 
(log ,) 3 -2.03) 


To the basic Hughes line, a form factor is then added. The purpose of this form factor 
is to help account for the differences in the laminar to turbulent transition seen by a 
flat plate and a fuller body. The attractiveness of Hugh's method 1s that it attempts to 
keep friction and eddy resistance linked. Physical observation of a bodv moving 
through the water would seem to make this reasonable. The method proposed by 
Hughes has not received the same acceptance within the Naval Architecture 
community as the ITTC 1957 line. The earliest criticism was the low value that 
equation (2-6) vields when compared to flat plate testing. Additionally. the use of a 
form factor is also questionable as outlined by Oosterveld’, quoted below: 


"It has now become clear that in many cases the wave 
pattern resistance is less than the residuary resistance estimated 


° Oosterveld. MWC. ed. "Report of the Resistance Committee”, 1 Sth International Towing Tank 
Conference, The Hague, 1978, p 2] 


by assuming a form factor independent of Froude number and 
Revnolds number. This will certainly be the case whenever bow 
wave breaking occurs. One possible reason for this 1s that the 
viscous resistance may varv significantly with Froude number 
(which effects the flow at the edge ot the boundary laver) as 
well as with Reynolds number, so that it 1s not well predicted on 
a constant form factor basis. However it 1s conceivable that this 
variation arises at least in part from an interaction of the 
boundarv laver and wake with the wave svstem, so that energy 
losses which originate with wave making, and are pnmanily 
Froude number dependent. may manitest themselves as an 
apparent increase in the viscous wake losses.” 


The tnctional lines discussed in this section are shown in Figure (2-2). 





75 
2008 ees ec eee aeee 
L065 R,-2) 
IN mencome #7" ine 2522-106, 1h, eC, 
007 Be wl, A 
, x 
YS awe ee BES LINE Spec O66 : 
XN 1LOG,, R,- 203! 
20061 Bs : 


-- ae ee me eee “> 





. —_- -4 ———— 


aos ; 


~ iTTC LINE ; 
etd ey ~ - 


DE ES Se a eee ee = eS 
.» 
= oe ATTC LINE «00004 
NPA | ! 

C009 ee = 

| aTTc Lin€ “Se Toa 

see, 

0002 = “~ 


Figure 2-2: Standard Skin Friction Lines’ 


This discussion of frictional resistance has beanng on the SLICE project in two 


° Pnnciples of Naval Architecture, op.cit. p. 298 
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respects. First there is the issue of turbulent flow. The general practice in model 
testing 1s to artificially stimulate turbulence. This is done so that the model and the 
full scale ship operate in the same regime. Without turbulence tripping, most models 
would operate at least partially in the transition regime shown in Figure (2-1), making 
data interpretation difficult since C, is highly variable in this regime. Hughes’ 
describes standard methods of stimulation (usually pins, wires or sand strips), 

methods of accounting for the resistance caused by these devices and the uncertainty 
involved in insuring that turbulence 1s actually achieved. The values of 8 at which the 
SLICE models and full scale ship operate are shown in Table (2-1). A comparison of 
these values to Figure (2-1) would appear to demand that some form of turbulence 


tripping be used in testing. 
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Full Scale 1/8th Scale 1/16th Scale 
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(Turbulent in (Transition in (Transition in 
Num ber Fig 2-1) Fig 2-1) Fig 2-1) 
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Table 2-1: Operating Reynolds Numbers for SLICE 
based on Strut and Pod Length 


Figure (2-1) however shows the laminar to turbulent transition for a flat plate. SLICE, 
with a nearly hemispherical nose shape and an overall L/D = 4.5 cannot be likened to 


a flat plate. The laminar to turbulent transition for a sphere is known to occur in the 


: Hughes.G and Allan J.. “Turbulence Stimulation on Ship Models. SNAME Transactions, 1951, 
p281-314 
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region NR =4.5 x 10°- 5 x 10° or two orders of magnitude lower than that of a flat 
plate. The actual transition for SLICE probably lies somewhere between these two 
extreme boundaries. When viewed in this light, turbulence stimulation is probably not 
worth the uncertainty involved in mathematically removing the resistance effects of the 
trip device. For the purposes of testing in this project. at the model scale. turbulent 
flow will be assumed even though turbulence stimulation is not provided. 

The second issue 1s to decide the most appropnate method to compute the 
frictional component of SLICE resistance. The ATTC line should be excluded since it 
is derived strictly from flat plate results. The Hughes line with form factor corrections 
may initially appear attractive since it incorporates calculated adjustments for form. 
The fullness of the SLICE body and the definite appearance of a breaking bow wave 
may cause the inaccurate accounting discussed by Oosterveld. The ITTC 1957 line is 
the best choice. [t 1s based on actual ship and model results. Even though there may 
be uncertainty in describing the details of laminar to turbulent transition, and the 
impact the fullness of a ship shape has on that transition. the ITTC 1957 line must be 
thought of as correctly including the phenomenon since it is based solely on test 
results. This is a fallout of the model to ship correlation method used to derive the 
line. It 1s the only line discussed that was developed in this fashion. Further work in 
this project will utilize the ITTC 1957 line as the basis for calculating frictional 


resistance. 


2.3 Wave Making Resistance 

Wave making resistance has two parts. In an obvious way it accounts for the 
energy expended creating the visible wave system seen as a bodv moves through the 
water. It also accounts for the energy expended overcoming drag due to the pressure 
distribution created by unseparated inviscid flow around a body. The literature is 
vague on the second part of this definition of wave making resistance. The pressure 
distribution around the body is directly related to the form of the body but this effect 
is not the commonly termed "form drag". "Form drag" is a Reynolds scaled viscous 


effect that is best thought of as the correlation between flat plate frictional resistance 
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prediction and actual ship friction results. The form part of wave making resistance 
should not be confused with eddy resistance either which accounts for the effects of 
separated flow. Since both parts of the definition of wave making resistance are 
pressure effects. this phenomenon will follow Froude scaling. 

One of the reasons that model testing is performed is because the ability to 
calculate wave making resistance remains an elusive goal. This section reviews the 
classical work that has been accomplished in this area. The objective is to gain some 
insights into the physical phenomenon and how it applies to the SLICE hullform. 
The foundation for a mathematical approach to wave resistance was set in place at 
nearly the same time that Froude was developing his model testing methods. J.H. 


Mitchell gave an analvtic solution for wave resistance of the form’: 
R = 08 [sec?0 | f foexpt(-&)sec?0(y-ixcos®)\dxdy /7d8 (=) 
wie Vy? 


Where p is the density of water. V is velocity and o is the potential function that 
represents the bodv of interest. The combined complexity of writing o to comply with 
the necessary boundary conditions at the free surface as well as with the relatively 
complex form of a monohull and the difficulty in actually solving this integral without 
the advantage of numerical computing techniques, were serious setbacks to the 
popularity of this approach. The most significant pursuit of Mitchell's work was 
undertaken by T. Havelock and W.C. Wigley. Wigley was the first to successfully 
decompose the waves generated by a hullform into components and thus gain an 
understanding of how thev are formed. This was first achieved with what 1s now 
called the Wigley hull. a wedge followed by parallel midbody followed by an identical 


wedge pointed opposite to the bow. Eventually Wigley refined the hull to a smooth 


* Newman. J.N., Marine Hydrodynamics. The MIT Press. Cambridge.1989. p.282. 
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form. He discovered that there are actually four waves produced by this shape. The 
first is produced by the leading edge of the hull, the next two are produced by the two 
shoulders formed by the intersection of the fore and aft wedge with the parallel 
midbody and the final one is produced by the trailing edge of the aft wedge. This 
decomposition held true for hullforms with continuous lines as well. Figure (2-3) is 
reproduced from Wigley's writings which graphically demonstrates this additive 
Sree 

This discussion is useful to the SLICE program because it gives insights to an 
approach at hullform optimization. Since the amplitude of the wave system generated 
is related to the form of the body, maintaining lines that are reasonably fine should 
help to minimize the hump of the wave making curve. As a approximation without 
actual calculations, the shape of the struts and pods of the SLICE form will be 
adjusted in an effort to find a reasonable compromise between frictional and wave 
making resistance. This is different from the LMSC approach which appears to have 
concentrated on the operating range of ¥ > 0.9 with effort focused on minimizing 
frictional resistance. 

The SLICE hullform can be decomposed into eight different bodies, four struts 
and four pods. Each of these component bodies has a form that is reasonably similar 
to the hullforms tested by Wigley. If each of these bodies can be thought of as 
producing four waves that combine into a single wave system for the body, then an 
approach to optimizing wave making resistance can be created by adjusting body 
spacing once form issues are decided. A three tier svstem of looking at this problem is 
proposed: 

1) First consider the interaction between a strut and its attached pod. This step 
can be thought of as creating an oversized bulbous bow for the strut. 

2) Next consider the effect of longitudinally aligned strut pod 


combinations. This will probably have a diminished effectiveness compared to step | 


. Wiglev.W C.S.. “A Comparison of Expernment and Calculated Wave-Profiles and Wave- 
Resistances for a form having Parabolic Waterlines", Proceedings of the Royal Society, Series A, 
London, 1934. pp 144-159. 
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Figure 2-3: Calculated Wave Profile Analysis 
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since the amplitude of the forward system will not equal that of the aft system once 
they interact. 

3) Finally consider the interactions of the transversely and diagonally aligned 
strut pod combinations. Since these are not in the principle direction of motion, it is 
expected that the effectiveness of these interactions will again be diminished from the 


previous step. 


2.4 Eddy Resistance 

As previously noted. the literature is noticeably vague in defining both form 
and eddv resistance. They are frequently used together but are not used 
interchangeably. Articles written about resistance prediction usually account for eddy 
resistance as an input to residuarv resistance. The research for this thesis uncovered no 
direct usage of form resistance 1n prediction schemes. As discussed in Chapter 2-2, 
this thesis considers form drag to be the increase in frictional resistance felt by a ship 
shape compared to an equivalent flat plate. The ITTC 1957 line directly accounts for 
form drag since it is based on experimental data from monohull forms. Other schemes 
such as the ATTC line and Hughes Line require some form compensation for accurate 
use since they are based on flat plate resistance. 

Eddy resistance on the other hand has physical origins. In this thesis it is taken 
to mean the result of separation effects from discontinuities in the hull form. These 
can result from actual breaks in the ships lines such as the intersection of the strut and 
pod at both the leading and trailing edge. intersection of appendages and the main 
body, or radii of curvature which are too hard for the flow to follow. As discussed in 
chapter 2-1, eddy resistance 1s a pressure effect and therefore should be scaled as part 
of the residuaryv resistance component. 

Following the classic approach, the minor swirling seen along the sides of a 
ship are associated with viscous effects of well behaved flow and are accounted for in 


the frictional component’’. Significant swirling that can be associated with flow 


Principles of Naval Architecture. op cit. p.312. 
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separation, such as that seen at the transom, describes the tvpe of effect accounted for 
by eddy resistance. There is some obvious room for interpretation of cause and effects 
in these definitions of form and eddy drag. This highlights the importance of observing 


physical behavior through model testing. 


2.5 Spray Resistance 

The formation of spray and the impact it has on total drag has historically been 
assumed insignificant to monohull testing. Simple observation of full and model scale 
monohulls operating at speed validates this assumption. Although spray 1s generated, 
the region of generation 1s highly localized and the quantity of spray generated is 
usually small. With the increased interest in advanced hullforms such as hydrofoils 
and planing craft. neglecting spray was recognized as a poor assumption but the body 
of available knowledge regarding the mechanism of sprav generation and the 
appropriate method of scaling spray results is relatively small. 

When spray becomes a major concern, the one equation. one unknown analogy 
used in section 2-] becomes one equation and two unknowns. The typical solution to 
this problem has been to test at verv high speeds (.¥ > 3.0) to eliminate wave making. 
This allows spray drag to be analyzed separately but is only a partial answer. A gap 
exists at the lower, more practical values of % In 1971, Chapman’ conducted a 
significant test program on strut forms with parabolic leading and trailing edges. He 
concluded the following: 

1. Spray drag is parually dependent on both strut form and thickness. Blunt 
leading edges produce more spray. 

2. The spray mass flow rate of spray depends on the thickness to chord ratio 
(t/c) and on .. Blunt bodies send the sprav up at higher angles causing greater 
wetting. 


3. Skin friction due to wetting 1s the primary source of spray drag. 


i Chapman. R.B. . Spray Drag of Surface Piercing Struts, Naval Undersea Research and 
Development Center Report AD 730710. 1971. 


4. The momentum in the spray sheet is sufficient to create a large pressure 
drag on anv object it strikes. 

Chapman's work was expanded by Zhu Bing-quan in 1986'*. A total of twelve 
struts were tested, all with parabolic trailing and leading edges. The objective was to 
correlate strut performance to the parameters c (chord length), ¢ ( maximum thickness) 


and x ( longitudinal position of the maximum thickness). The relationship derived is: 
_ 12,2 
R=pV2x2C, (2-8) 
and C, 1s given as: 


C,=0.068- -0.555(=)?+0.696(=) 0.0832 F (2-9) 


The research also showed that the spray drag was minimized when ¢.x = 0.313. 

This has applications to the SLICE analysis. The strut form selected by LMSC 
is contrary to that recommended by the published literature for minimizing spray drag. 
The apparent design goal of minimizing frictional resistance fails to consider spray 
effects. In deciding the extent to which the sprav effects will scale from model size to 
full scale. the dependence of (2-9) on ¥ implies that sprav drag will follow Froude 
scaling. The conclusion reached by Chapman and Zhu Bing-quan that the initiation of 
spray 1s highly form dependent also supports this assertion. 

The following is recommended when dealing with spray effects. Choose strut 
shapes that are known to minimize spray generation. Unless there is a significant 
structural or interior arrangements reason to do otherwise, struts should have parabolic 


leading and trailing edges and / x should be reasonably close to 0.313. For the ATD, 


? Zhu Bing-quan, Ge Wei-zhen. “Experimental Study on Spray Performance of Surface Piercing 
Struts, Journal of Hydrodynamics, ser 2 no 4. 1990, pp 91-98. 
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this will help minimize spray effects and allow the research efforts to focus on the 
area of principle interest, the wave making resistance. Anv spray that is observed in 
testing should be treated as part of the residuary resistance. If the tentative conclusion 
that spray follows Froude scaling is wrong, then this process will over predict the total 
resistance. Considering the uniqueness of the SLICE hullform, a conservative error of 
this nature 1s well justified for the first full scale ship. Design work in this project will 


follow these conclusions. 


2.6 Appendage Resistance 

The effect of appendages on total drag has historically been assumed small in 
monohull testing. This assumption 1s based on the fact that the overall size of the 
appendages 1s usually one or two orders of magnitude less than the main hull. For 
comparative testing of hulls. 1t 1s reasonable to ignore the appendages. For obtaining 
exact resistance predictions, appendages are usually included in the final testing and 
stripping tests have been devised to help account for their effects. In stripping tests, 
the model is towed first with all the appendages on. In subsequent tests, the 
appendages are removed one at a time. In this manner the impact of each appendage 
can be assessed. This approach assumes that the interaction effects of the appendages 
with each other and with the main hull are negligible. 

For SLICE, appendages occupy a larger percentage of total volume compared 
to monohulls. The flow around SLICE is more complex and therefore interactions may 
be more pronounced. A full stripping test 1s warranted in the case of SLICE. Since 
this project 1s focused a comparison of alternate hull forms, appendages will be 


excluded, leaving the stripping test to future work. 


2.7 The SLICE Data Scaling Method 

Based on the above discussion, the traditional method of data scaling outlined 
in section 2-1 is recommended for the SLICE project and will be used in this research. 
Calculating the frictional resistance using the ITTC 1957 line will yield a moderately 


conservative answer. Although the ITTC 1957 line was denved from monohull tests, it 
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1s the only method that naturally accounts for the Froude Number dependence of 
boundary laver development. All other phenomena (wave making, eddy formation and 
spray development) will be treated as part of residuary resistance and scaled 
accordingly. 

It is important to recognize is that this technique of data scaling is considered 
reliable due to the enormous volume of model testing that has been performed on 
numerous series of monohull models. The validation of this process for SLICE will 
only come when the full scale ship 1s built and tested. This is the real value of the 
ATD. Building and properly testing the full scale ship will significantly contribute to 
the Naval Architecture community's understanding of the robustness of current testing 


and scaling techniques. 


2.8 A Numerical Approach 

Model testing is a time consuming, detail intensive process that requires 
significant capital expenditure. For a new hull form. some degree of model testing is 
mandatory in the early stages of design so observations of phvsical behavior can be 
made. This is critical to validating assumptions, but from a cost savings perspective, 
the ability to numerically model the ship's behavior is attractive. Small changes to the 
hull can be made and tested without the expense of building a new model. If the 
numerical model can be kept mathematically simple, comparative results can be 
obtained quickly and cheaply. There must be a degree of confidence however that the 
computer model fully represents actual resistance. 

The phvsical phenomena of interest are friction, wave. spray and eddy 
resistance. Friction and spray resistance can be calculated using (2-5), (2-8) and (2-9). 
If there is no significant separation or vortex phenomena, it is reasonable to neglect 
eddy resistance (only observation can answer this) and a prediction for full scale 
resistance can be obtained if a method of solving (2-7) can be found. This is 
significant to the SLICE discussion since it may provide a simple and fundamental 


tool for predicting and minimizing the wave making hump magnitude. 


A solution to Mitchell's integral was outlined by Lunde” in 1951. The process 
was further applied to a ship called the Wave Cancellation Multi-hull Ship Concept by 
Wilson and Hsu’ They express the total wave making resistance as the summation of 


body and interaction terms as shown: 


Bis 
NN N 


2 
: 2-10 
R, =16n px’ [LEP P+ L0,.Q)sec"@d8 (2-10) 
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where N represents the number of bodies the overall hull can be divided into and: 
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'* Lunde. JK. “On the Linearized Theory of Wave Resistance for Displacement Ships in Steady 
and Accelerated Motion”. SNAME Transacnons, Vol 59, 19581. pp. 25-76. 


'’ Wilson. MB. and Hsu. C C., "Wave Cancellation Multhull Ship Concept", Presented at the 
Intersociety High Performance Marine Vehicle Conference. Arlington, Va.. 24-27 June 1992. 
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In these expressions, the like product terms are the body terms and the cross product 
terms are the interaction terms. 

The underlying assumptions to Lunde's approach are : 

1. The fluid is incompressible. trrotational and inviscid. 

2. The wave height is small compared to wave length. 

3. The ship motion can be reasonably represented by a fine distribution of 
sources and sinks along the centerline of the bodv. 

Prior observation of SLICE models in the towing tank indicates that conditions 
2 and 3 mav not be satisfied ideally. Each of these phenomenon must be watched 
closely in subsequent tests and considered in the interpretation of results. Havelock’? 
demonstrated good correlation using this technique with hull forms having L/B as 
small as 10.6. Data presented by Newman'® indicates that there is better correlation 
between model test data and thin ship theory at higher values of ¥. SLICE is a full 
form with L.D = 4.5-5.0 but operates at ¥% = 0.9-1.2. Although SLICE is not an exact 


fit for thin ship theory, the approach mav be reasonable. 


'S Havelock, TH., "Wave Resistance Theory and its Application to Ship Problems", SNAME 
Transacnons, 195], pp. 13-24 
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Chapter 3 Model Testing 


As part of this research. three models were designed. built and tested. Appendix A 
provides detailed lessons learned from this testing experience which may be useful to 


other students attempting to undertake similar model testing projects. 


3.1 Model Test Objectives 
There were three objectives underlving the model testing phase of this project: 

1. Expand the quantity of available model scale resistance data for the 
baseline SLICE configuration. 

2. Visually observe model performance in sufficient detail to 
understand the physical behavior of this hullform. As discussed 1n chapter 2. this 1s 
important for justifving the scaling technique selected. 

3. Test operationally equivalent hullforms to help assess baseline 
SLICE advantages and disadvantages. 

The models were tested at the United States Naval Academy's 380' tow tank 


facility during the week of 6-10 March 1995. 


3.2 Model Descriptions 

The three models tested as part of this research were designed to meet the 
same operational constraints. This was done to ensure that comparisons made at the 
end ot the testing could be made on the basis of similar ships. All models when scaled 
to full scale have a displacement of 177 Lton. a pod diameter of 7-8 feet and a strut 
beam of 3.25 feet. These dimensions are the working dimensions of the ATD vessel, 
driven by requirements other than hvdrodvnamic performance. The pod diameter 1s 
controlled by engine size requirements and the strut beam 1s controlled by access 
trunks into the pods. With these requirements imposed on all three models, they are 
considered equivalent hulls. 


The three models tested were designated M-1. M-2 and M-3. Detailed tables 
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of hull offsets and design calculations for each model are contained in Appendices B. 


C and D respectively. Individual descriptions of each model follow. 


MODEL M-| BASELINE: 

This model utilizes LMSC's offsets as tested in December 1994 with the 
exception that the fillets at the strut/pod interface and the haunch structures are 
excluded. The fillets were excluded to help simplify model building. This modification 
is not considered detrimental to the testing outcome. The combination of the pod, strut 
and fillet in the LMSC model creates a effective L/D of 3.9, the omission of the fillets 
increases this ratio to a more favorable 4.2. This modification should result in slightly 
lower resistance values. The bow wedges installed in the December LMSC model are 
replicated in this model. The purpose of this model 1s twofold. First. it provides a 
means of correlating the data obtained in these tests with the data obtained by LMSC 
in December 1994. This 1s considered important since the tests were conducted in two 
different facilities. under different conditions with different scale models. The ability 
to correlate data between Model M-1 and the LMSC model will be valuable in 
demonstrating repeatability of results. The second use of Model M-1 1s to expand the 
available range of data for the baseline SLICE. Model tests were conducted at model 
equivalent speeds through 40 knots to ensure that the entire operating range was 
properly represented by data. 

Table 3-1 provides the principle dimensions of Model M-1. Profile and top 
view sketches of the forward and aft strut/pod combinations are in Figure 3-1. In this 
drawing the alignment of the struts with the corresponding pods 1s correctly 
represented. The spacing between the forward and aft combination 1s not preserved 
however in the interest of fitting the drawing onto one page. This model was 
constructed with 29 Ibf of reserve buovancy. The struts are manufactured from balsa 


wood with a fiberglass shell. the pods are molded fiberglass. 
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Table 3-1: Model M-1 Baseline 
Principle Dimensions 
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Figure 3-1: \Mlodel \l-1 Forward and Aft Strut/Pod 
(Note: In this drawing the fwd-aft pod spacing 1s not preserved) 
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MODEL M-2 RATIONAL DESIGN: 


The purpose of this model is to explore the effects on total resistance of using 
finer offsets for the struts and pods . To avoid changing too many variables in a single 
test, this model maintains all the gross parameters of the Model M-1 Baseline (ie: 
Strut/Pod alignment, fore/aft spacing and centerline to centerline spacing). The 
expectation for this design 1s that by comparing this model's results with that of M-1, 
some conclusions can be drawn about the relative magnitude of both friction and wave 
making resistance through the range of speeds tested. The anticipated result from this 
model is that the amplitude of the wave making hump will be diminished. 

The offsets for the pods are adapted from Series 58 data'’. This series was 
selected since it represents a concerted effort to optimize submarine hull forms. The 
model selected from the series was Model 4155. This model was selected based on its 
length to diameter (L/D) ratio of 5 which 1s similar but finer than the baseline SLICE 
L/D of 4.2. It should be noted that within the Senes 58 results, Model 4155 falls 
outside the range of optimum L D which was found to be between 6 and 8. Model 
4155 was selected despite this drawback in order to satisfy the operational requirement 
of creating a full scale 177 Lton vessel with pod diameters of 7-8 feet. These critena 
could not be satisfied with a senes 58 hull with a L/D of 6-8. 

A trait of the Series 58 hullforms is that they have no parallel midbody. 
Although hvdrodyvnamically this 1s superior. it would significantly complicate the 
building of the full scale ship. The Model M-2 therefore incorporates a modification to 
Model 4155 by replacing some of the midships shape with parallel midbody. 

The strut offsets utilize a symmetric parabolic leading and trailing edge. The 
general parameters of this strut are outlined in Chapter 2-5 but are also constrained by 
the design requirement that the full scale maximum beam must be at least 3.25 feet. 
The calculations used to develop the strut and pod offsets and a table of the actual 


'’ Gertler. M. "Resistance Experiments on a Systematic Series of Streamlined Bodies of 


Revolution for Application to the Design of High Speed Submarines", David Taylor Model Basin Report 
C-297. unclassified and approved for public release 18 Dec 1972, April 1950. 
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offsets used are contained in Appendix C. 

Model M-2 represents a partial optimization. The offsets for the struts and pods 
are based on research conducted in optimizing these tvpes of shapes. The optimized 
values are then adjusted to conform to design requirements. Model M-2 make no 
attempt to optimize the placement of the strut pod combinations relative to one 
another. Table 3-2 provides the principle dimensions of Model M-2. Profile and top 
view sketches of a the forward and aft strut/pod combinations are at Figure 3-2. This 
model was constructed with 30 Ibf of reserve buovancy. The struts are manufactured 


from balsa wood with a fiberglass shell, the pods are molded fiberglass. 


MODEL M-s OPTIMIZED HIGH SPEED SWATH: 


The purpose of this model 1s to compare SLICE performance to a more 
conventional SWATH. The offsets used for the high speed SWATH were reported by 
McGregor’®. The demihull diameter was forced to meet the 7-8 foot full scale 
requirement. the strut maximum beam was forced to meet the 3.25 foot full scale 
requirement and the overall displacement was forced to meet the 177 Lton full scale 
requirement established earlier. In order to meet these requirements, the demihull L/D 
was established at 12.86. This corresponds to a full scale demihull length of 90 feet, 
about 15 feet shorter than the Length Overall (LOA) of the proposed ATD SLICE. 
The consequence of this is that although Model M-3 is hydrodynainically equivalent to 
Model M-1. depending on the extent to which the box structure 1s cantilevered forward 
and aft of the struts, the SWATH will end up with as much as 600 square feet less 
arrangable deck space. This 1s an 11% decrease from the SLICE ATD vessel. 

A tandem strut arrangement was selected in order to minimize the displacement 
occupied by the struts. This was necessary in order to maximize the demthull 


displacement and maintain a reasonable demihull length. As with model M-2, the 


‘McGregor. R.C. and Chun, H.H., “On the Potennal of SWATH Ships for Very High Speed 
Operation”, Fast 91, pp. 491-506 
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strut offsets are based on the guidance provided in Chapter 2-5. In his SWATH 
design, McGregor provided guidance for the optimum placement of the struts in 
relation to the pods. This placement was utilized in Model M-3. The calculations used 
to develop the strut and pod offsets and a table of offsets are contained in Appendix 
D. Table 3-3 provides the principle dimensions of Model M-3. A profile and top view 
sketch of a the demihull 1s at Figure 3-3. The struts of this model are manufactured 
from balsa wood with a fiberglass shell, the demihulls are hollowed PVC rod, shaped 


to the proper offsets. There 1s 10 lbf of reserve buovancy. 


The reason the three models were not all built to the same scale factor stems 
from the actual building process. There 1s a slight vanation in the full scale diameter 
of the three underwater hull diameters. In full scale. M-1 is 8 feet. M-2 is 7.74 feet 
and M-3 1s 7 feet. PVC pipe was used to mold parallel midbody sections. This 
material 1s relatively inexpensive and comes in a_ standard pipe outer diameter 
dimension of 5.5625 inches. This dimension became the dnver for determining the 
model scale factor. Since the three models have different full scale diameters, the scale 


factors varv accordingly. 
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Table 3-2: Model M-2 Modified Baseline 
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Figure 3-2: Model M-2 Forward and Aft Strut/Pod 
(Note. In this drawing the fwd-aft pod spacing is not preserved) 
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Demthull Length 


Strut Length 19.10 inches 


Strut Max Beam 2.58 inches 





Table 3-3: \Miodel 1-3 McGregor's High Speed SWATH 
Principle Dimensions 





Figure 3-3: Model 1-3 Demihull and Strut 


+] 


3.3. Test Description 

A copy of the test plan developed for the Naval Academy Tow Tank tests 1s 
included as Appendix E. This enclosure explains the planned tests in detail. 

Two basic tests were planned for the models. The first measures resistance and 
lift with the models fixed in pitch and heave. This test replicates the test configuration 
used by LMSC thereby allowing for correlation with the LMSC data. Data collection 
runs consisted of the model held at its design draft with no trim. At the slow end of 
the speed spectrum, one speed change was allowed during data collection. At higher 
speeds only one speed was used per test run. Data was collected at full scale 
equivalent speeds from 3-40 knots with the following model configurations: 

1) Model M-] Baseline in the LMSC configuration. 

2) Model M-1! Baseline with the forward and aft struts in line. 

3) Model M-1] Baseline in the LMSC configuration at varving drafts. 
4) Model M-2 Modified Baseline in the LMSC configuration. 

5) Model M-3 McGregors SWATH. 


The second test, quasi free in heave, was not conducted. The intent of this test 
was to replicate the effects of control surfaces through a counterweight system with 
the hope that information could be gained regarding the necessary sensitivity of the 
control system. Conceptually, the counterweight would provide the same lift as the 
control surfaces without added model building complexity. The model would then be 
free to respond to small perturbations and its stability qualitatively assessed. This test 
was considered a second priority to the project and time in the test facility ran out 
before it could be accomplished. It is considered a worthwhile area of further study. 

The test plan indicates that a fourth model. the Midfoil, was tested. This model 
was built and tested at the request of ONR but is not considered part of this research. 


It is not reported in this thesis. 
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3.4. Physical Observations 

The importance of observing the behavior of the models can not be overstated. 
A videotape of the model tests is included as part of this thesis. Shown in the tape is a 
port beam perspective, both above and underwater as well as a stern view of all the 
models tested. The video is organized by speed increment, showing each model first at 
the peak of the wave making hump (10 and 12 knot full scale equivalent) then in the 
post hump valley (16 knot full scale equivalent) and finally at three high speed 
increments (25, 30 and 35 knot full scale equivalents). The following qualitative 


observations are made regarding the performance of the models. 


PROPELLER SUBMERGENGE 

The ATD vessel will have propellers mounted on the aft end of the forward 
pods. After the August LMSC model tests, there was significant concern that the wave 
system created in the hump was large enough to expose the propeller. This evaluation 
was made bv observing the wake from above the free surface but was difficult to 
assess properly due to spray and distortion. With M-1, the underwater camera used in 
this project reveals that there is better submergence at the propeller than apparent 


when viewed from above. 


HUMP SPEED PERF@RMANGE 

ft was expected that M-2 with its finer lines would exhibit less wave making 
than the very bluff M-1! in the hump speed range. In fact, the surprising result was 
that the wakes produced were nearly identical. The wake of M-3 at the same speed 
was much smaller than either M-1 or M-2. The implications of this observation are 
unclear without examining actual data. The similarity between the M-1 and M-2 wakes 
would suggest that the form of the bodies, both strut and pod, is not a controlling 
factor in hump speed resistance. The vast difference between M-1/M-2 and M-3 
suggests that the underwater body plays a role in wave cancellation. 

The correct conclusion may simply be that this is a demonstration of the flow 


differences of bluff bodies (M-1 L/D = 4.2, M-2 L/D = 5.0) and more slender bodies 
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(M-3 L/D = 12.9) or there may be some more complicated mechanism of wave 


cancellation that needs to be explored. 


HIGH SPEED PERFORMANCE 

In the full scale range of 25-40 knots, there is a significant difference in the 
performance of M-1] and M-2. At these speeds, M-1 creates a verv large vortex at the 
aft intersection of the strut and the pod. Above the surface, this 1s manifest as a 
rooster tail. As the speed increases. the inception point of this vortex travels further aft 
along the pod tail and at the full scale equivalent of 40 knots almost reaches the end 
of the pod. This effect 1s absent in M-2 and M-3. 

This effect does not appear to be caused by the relatively small scale of the 
models. Reviewing footage of the LMSC tests, the same phenomenon 1s apparent. It 1s 
more pronounced in the August tests than the December tests, but it is difficult to tell 
if that can be attributed to the design changes LMSC made. There is not enough 
quality video footage of the December tests to draw this conclusion. It 1s possible that 
the August and December tests did not allow a full steadv state condition to be 
established. In both the August and December tests. speed increments were maintained 
on average for 10 seconds measured from the start of a speed change to the start of 
the subsequent speed change. As observed in the March tests, this amount of time 1s 
barely adequate to establish steady state conditions. When accelerated from rest. the 
March data shows that forces did not stabilize until 8-10 seconds after the steady state 
speed was achieved. The LMSC testing may have stabilized sooner since the 
acceleration was only between speed increments and not from rest. Nonetheless, the 
data collection period was extremely short. The inception speed of the rooster tail 
appears to be comparable in both the LMSC tests and this project's tests. This 
indicates that the effect follows Froude scaling. 

After the August LMSC tests. much attention was focused on scaling spray 
effects. This was in response to the extreme wake created at high speed. LMSC's 
position was that spray effects could be ignored in full scale predictions and designed 


out using spray rails. The underwater observation of M-1 demonstrates that much of 
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what was called spray in August was actually produced by an entirely different 
mechanism. The vortex that creates the rooster tail is a separation effect created as the 
flow unsuccessfully attempts to conform to the hullform. When M-1, M-2 and M-3 are 
compared at the same speed. it becomes obvious that only a very little of the splashing 
seen above the surface in M-1 can be attributed to spray generated by the struts. The 
majority of the splashing effect is more accurately accounted for as a by-product of 


eddy resistance. 


WAKE CHARACTERISTICS 

The wake produced by the SLICE hullform is compact and has a greater 
amplitude than that of the SWATH. The SWATH produces a relatively flat but 
significantly longer wake. From an energy in the water standpoint, comparison 


becomes difficult since the entire volume of the wake must be considered. 


3.5 Test Data 

Raw data collected during the testing 1s presented in Appendix F. Test | 
results (fixed in pitch and heave measurements of drag and lift) are presented 
graphically in. Drag (positive x) 1s taken to be opposite the direction of forward 
motion. Sinkage (positive z) 1s down from the free surface. Several speed settings 
were repeated during the testing. Repeated speeds are recorded in these figures and are 
designated Dual Point (DP). No direct comparison of performance using Figures 
(3-4) - (3-6)should be done since each model is built to a slightly different scale 
factor . 

Figures (3-7) and (3-8) present model scale data at the same scale factor. This 
required scaling M-2 and M-3 data to the M-1 scale. Data scaling was accomplished 
using the process outlined in Chapter 2-7. Scaling between models was accomplished 


using the composite scale factors shown: 
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1,,= Anodes = 12-19 <9 975 (3-2) 
Meee wee O20 


The calculations used 1n this process are in Appendix G. 

Figure (3-9) shows the data collected for Model M-1 tested in both the 
Baseline configuration and an in-line configuration of the forward and aft strut/pod 
sroups. The results of testing Model M-! at varving depths was inconclusive. 


Insufficient data was collected to develop a trend. 
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Figure 3-4: Model M-1 Fixed in Pitch and Heave 
Unscaled Data 
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Figure 3-5: Model M-2 Fixed in Pitch and Heave 
Unscaled Data 
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Figure 3-6: Model M-3 Fixed in Pitch and Heave 


Unscaled Data 
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Figure 3-7; Models M-1/M-2 and M-3 Fixed in Pitch and Heave 
Data Presented in 17.26:1 Scale 
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Figure 3-8: Models M-1/M-2 and M-3 Fixed in Pitch and Heave 
Data Presented in 17.26:1 Scale 
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Figure 3-9: Models M-1 Fixed in Pitch and Heave 
Tested in the Baseline and In-Line Configurations 
Data Presented in 17.26:1 Scale 
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Several comments can be made about this data. 


REPEATABILITY OF RESUES 

A large number of data points were repeated through the course of the testing. 
Without exception, the difference in recorded data is statistically insignificant. Most 
repeated data was recorded in successive runs down the tank. A notable exception to 
this are runs 70 and 71] in Appendix F. In this case. testing of Model M-1 was 
completed. Model M-2 and Model M-3 were tested, then Model M-] was remounted 
and realigned on the carriage and five data points collected. There was excellent 


correlation between these points and data collected two davs earlier. 


HUMP SPEED MODEL PERFORMANCE 
The resistance values for M-] and M-2 through the hump, Figure 3-7, support 
the observation made earlier that the wakes are very similar. The similarity in 
resistance values was an unanticipated result. M-2 represents the extent to which L/D 
can be extended and still retain the features of SLICE. Since this vanation was 
insufficient to impact the amplitude of the hump. the only other alternative for 
diminishing the hump magnitude is altering the relative position of struts and pods. 
The fact that the hump for M-1, M-2 and M-3 all occur at essentially the same 
speed indicates that the Strut Length is the appropriate length to charactenze SLICE. 
This conclusion is reached by considering M-3 in contrast to M-] and M-2. The hump 
is expected to occur at .%= 0.45. For the demihull of M-3 this corresponds to 
V= 6.25 fps. for the struts it corresponds to V= 3.2 fps. The hump actually occurs at 
V=4,5, 
demonstrating dominance by the struts. The pods of M-1 and M-2 are less than half 
the length of M-3's demihulls and vet the hump for all three is at the same speed. The 
struts for all three models were essentially the same length. The underwater bodies 
may be contributing to wave making through cancellation effects but the position of 


the hump 1s driven by the surface piercing bodies. 
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HIGH SPEED MODEL PERFORMANCE 

In the range 6.5-10 fps (full scale equivalent: 16-25 knots), M-1 total resistance 
IS up to 20% greater than that of M-2 even though M-1 has less wetted surface and 
therefore less frictional resistance than M-2. This is the first indication that model 
form must be traded off with fnction to optimize post hump resistance. 

For speeds greater than 10 fps. the rapid increase in M-] resistance can be 
directly correlated with the appearance of the trailing edge vortex noted in the physical 


observations. 


M-i LMSC COMPARISON 

One test objective for M-1] was to compare M-1 resuits with the December 
LMSC model. Figure 3-10 presents the M-1 data as collected as well as a sampie of 
LMSC data from both August and December scaled to 17.26:1. Model M-2 results are 
aiso presented scaled accordingly. The August LMSC curve represents LMSC data 
runs 106-i21, The December LMSC curve represents LMSC data runs 86-97. These 
modeis were tested with appendages and haunch structures in place. 

M-1 correlation with the December tests is good below 8.25 fps (full scale: 20 
knots). becomes marginal from 8.25-10 fps (full scale:20-25 knots), and does not exist 
above this value. M-1 correiation is more satisfactory with the August test results. 
The tack of high speed correlation requires comment. 

Of concern is the vortex created at 10 fps (full scale: 25 knots). The 
phenomenon is documented as existing in the August and December tests as well as 
the tests conducted on M-1. The question is deciding if the lack of correlation 
indicates that effect does not scale, or if there is some other mechanism that is causing 
the difference. The possibility that the LMSC data may have been obtained before 
steady state was achieved has already been discussed. Although this is a fault of the 
LMSC test program, it 1s difficult to attmbute the entire difference in test results on 


this snortcoming alone. 


(sh 
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Figure 3-10: Models M-1 Performance Compared With LMSC Results 
Data Presented in 17.26:1 Scale 
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In the tests of M-1. there 1s direct impingement of the rooster tail on the model 
frame. During the testing. this was considered acceptable for two reasons: 1) 
Impingement also occurred in the LMSC model tests and will occur in the full scale 
ship depending on the construction of the box structure. 2) During the testing, the 
impingement was considered a minor effect since the profile of the frame structure 
was small (60 in total). Reviewing video footage of these tests. the direction of 
rooster tail flow was perpendicular to the frame, thereby maximizing the force of 
impingement. In the LMSC tests impingement was not as direct. Appendix H 
provides a calculation that demonstrates that impingement could account for up to 
11 pounds of the measured force at 15 fps. 

There are enough differences in test configurations to preclude a definitive 
explanation for the lack of correlation. The rooster tail impingement may contribute to 
the lack of correlation. The conclusion drawn at this point is that the trailing vortex 1s 
real and will scale but may not directly contmbute to overall resistance to the extent 
indicated by the Model M-1! results in Figure (3-7). Failure to achieve steady state in 
the LMSC test may also be a significant factor in this difference. Re-testing M-1] with 
the frame isolated will help eliminate this ambiguity but ultimately full scale test 


results are required. 


3.6 Full Scale Predictions and Discussion 

With the understanding that the M-1 results may over predict full scale 
resistance, figures (3-11) and (3-12) present full scale data for all three models. Figure 
(3-135) provides power requirements based on these values. Data scaling was 
accomplished using the process outlined in Chapter 2-7. The calculations used in this 
process are at Appendix G. Based on figure (3-11), the high speed SWATH offers the 
best resistance performance through the entire speed range. 

Figures (3-14)- (3-19) present M-1. M-2 and M-3 data in model and full scale, 
decomposed into the frictional. residuary and correlation factor components. In the 
post hump region of the resistance curve. residuary (ie wave making) resistance 


accounts for at least 60% of the total resistance value. This further supports the 


56 


concept that the hull form must be designed as a compromise between frictional and 


wave making resistance, not a single component optimization. 
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Figure 3-11: Models M-1/M-2 and M-3 Fixed in Pitch and Heave 
Full Scale Resistance Data 
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Figure 3-12: Models M-1/M-2 and M-3 Fixed in Pitch and Heave 
Full Scale Resistance Data 
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Figure 3-13: Models MI-1/M-2 and M-3 Fixed in Pitch and Heave 
Powering Data 
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Figure 3-14: Model M-1 Resistance Components as a Percent of Total Resistance 
Derived from Unscaled Test Data 
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Figure 3-15: Model M-1 Resistance Components as a Percent of Total Resistance 
Based on Full Scale Results 
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Figure 3-16: Model M-2 Resistance Components as a Percent of Total Resistance 
Derived from Unscaled Test Data 
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Figure 3-17: Model M-2 Resistance Components as a Percent of Total Resistance 
Based on Full Scale Results 
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Figure 3-18: Model M-3 Resistance Components as a Percent of Total Resistance 
Derived from Unscaled Test Data 
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Figure 3-19: Model M-3 Resistance Components as a Percent of Total Resistance 
Based on Full Scale Results 
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Chapter 4 Numerical Solution 


4.1 Objective of a Numerical Approach 

The ability to predict a ship's resistance using numerical methods has many 
advantages. If an appropnate method can be found, the expense and time of model 
testing mav be reduced. Current state of the science 1s approaching realistic treatment 
of the free surface (1e: non lineanzed theory) but these methods are still too expensive 
to use as an iterative design tool. Simpler less expensive methods may not be accurate 
enough to exactly predict resistance but can be extremely useful for comparative 
analvsis of small changes to the same hullform. These methods require tow tank 
model testing as verification to ensure that all phvsical phenomena are properly 
represented. 

This thesis attempts to validate the SLICE applicability of a numerical method 
solution of Mitchell’s integral as an approximation for wave resistance. As cited in 
Chapter 2. the approach. first outlined bv Lunde. was used by Wilson in the Wave 
Cancellation Multihull Ship concept. It is not expected that this thin ship 
approximation will provide exact resistance values. The goal of this chapter is to 
develop an understanding of the body interactions in the SLICE hullform to assist 
resistance optimization through body placement. Of specific interest is reducing the 
magnitude of the wave making hump. The thin ship approach 1s useful since it allows 
wave making resistance to be segregated into all component body and interaction 
terms. This is important while studying the effect of small geometry changes on total 
resistance. 

Model M-2 is selected as the subject for this study. Physical observation during 
the model tests demonstrate that this model has a reasonable chance of satisfying thin 


ship flow requirements. 
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4.2 Solution Technique 


Restated from Chapter 2.8. the equation to be solved is: 


s NN NWN 
R,,, =l6n px | (LLP P +LLQ.0,)sec*ad6 
W 11 Tie 


where: 


TL 


p= { a(x,z)cos[K(xcos® +ysin®)sec6 Je =" Sard: 
00 


TL 
Q= / if 0 (1,2)sin[k(xcos® +ysinB)sec*B]e =** dud 
00 
o=-—— (fx2)) 
2m Ox 


3 
y- 


(4-1) 


(4-la) 


(4-1b) 


(4-1c) 


(4-1d) 


The coordinate system that will be used has its ongin in the undisturbed free surface 


with +X in the direction of ship motion. +Z vertically up from the free surface and 


+Y in accordance with the nght hand rule. 


Br ELOPNIE VW? OF THE BODY AND POTENTIAL FUNCTIONS 


The first requirement is to wnite a function that descnbes the strut and pod 


forms. This function is wnitten with the bodies at the global ongin. Later. these bodies 


will be translated to their actual relative positions through the X, Y and Z terms in 
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Equation (4-la) and (4-1b). 


For M-2. the strut form was developed using an identical parabolic nose and 


tail shape with parallel midbody in between. The equations used to describe the strut 


shapes are: 


Nose Offers == (1 -(“2)""™] 


MidbodyOffsets= 


to | & 


Tailogfsers == +{1 (=) 


where for this specific strut. 


B (the maximum beam) = 0.195 ft 
L, (nose length) = L, (tail length) = 0.51 ft 
Lp (parallel midbody length) = 0.42 ft 


= 795 
te = 22S 


(4-2) 


(4-3) 


(4-4) 


To reduce the complexity of the strut body potential function as applied in (4-1a) and 
(4-1b), the offsets calculated by (4-2), (4-3) and (4-4) were curve fit to a single high 


order polvnomial given bv: 


Strut =0.4893 *x°+2.1126 +x°+3.1328 «x* + 1.7248 +x°-0.0786 +x? -0.3890 +x +0.0008 
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(4-5) 


The source distnbution representing the strut is wnitten as: 


=- 2 (2,936 +x5+ 10.563 ¢x*+12.531 +x? +5.174 #x?-0.1572 +x-0.389) 
a It 


9 se 


The pod shapes are created in a similar fashion. Now however. the nose is 


represented by an ellipse and the tail bv a parabola. The pod offsets are given by the 


following: 


Lo=y 1, SS 
PodNose= = +(1-(——) ow (4-7) 


A 


PodMidbody=— (4-8) 
-f -j n 
PodTail=— +[1-(——* 2) - (4-9) 


t 


where for this specific pod. 


D (the maximum diameter) = 0.4635 ft 
L, (nose length) = 0.67 ft 

I Gail length) = 1.1 ft 

Lop (parallel midbody length) = 0.60 ft 
Meee 2.00) 

ey a 


Again a single high order polynomial was generated to represent the pod and is given 


by: 
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Pod=-0.04 +x ?-0.32 +x ®-1,16+x°-2.37 +x*-2.96 +x3-2.27 +x7-0,99 +x+0.038 


The source distribution representing the pod is written as: 


Cg *(-0.253 +x°-1,93 +x°-5.82 +x*-9,49 +x3-8.88 +x7-4.54 +x-0.99 
Dar 


The pods are treated as a line source in (4-1a) and (4-1b) bv fixing the z term 
in the exponential as a constant equal to the strut depth plus the maximum radius of 
the pod. This follows the example provided by Chapman’? for dealing with submerged 


bodies ot revolution. 


PAIRING OF TERATS 

Struts and pods are numbered in the traditional naval manner. The scheme is 
presented in Figure (4-1). There are a total of 36 possible combinations of the eight 
bodies that form the SLICE hull. Eight combinations are “like” or bodv terms and 28 
combinations are different or "cross product" terms. These terms are presented below 
grouped bv relative position in the hull. The numbers indicate which bodies are 
interacting in the term calculated: 

1) Bodv Terms: 
1-1. 2-2, 3-3, 4-4, 5-5, 6-6, 7-7. 8-8 


2) Same Corner Cross Terms: 


vs Chapman. R.B., Hydrodynamic Drag of Semisubmerged Ships, American Society of Mechanical Engineers. 
paper number ~2-Wa, 1972. pp 879-884. 


aa 


1-5, 2-6. 3-7. 4-8 
3) Longitudinal Cross Terms: 

1-3, 1-7. 5-3. 5-7, 2-4. 2-8. 6-4. 6-8 
+) Diagonal Cross Terms: 

let. 1-8. 5-4. 5-8. 2-3. 2-7. 6-3. 6-7 
5) Transverse Cross Terms: 

1-2. 1-6. 5-2. 5-6. 3-4. 3-8. 7-4. 7-8 


This physical grouping 1s retained through the calculations to simplifv companisons. 
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Figure +1: SLICE Body Numbering Scheme 


CODING 

MATLAB was used for the required computations. This choice was made to 
utilize the installed graphics capabilities of the program. Plotting results as 
computations were made was useful both in validating geometry assumptions and 
debugging. A logic flow chart of the overall computation scheme is at Appendix I. 
The routines used are also presented in this appendix. The pod and strut geometnes 
were developed in the files PODGEO.m and STRUTGEO.m, The main wave making 
resistance routine is contained in the file Q.m. The solution to equations (4-1), (4-1a) 
and (4-1b) are solved in called subroutines for each of the thirty six paired terms in 
files named for the bodies that are paired. Samples of a strut bodv, pod body and 
strut/pod cross term are contained in the files ONEONE.m. SLXSIX.m and 
ONESIX.m. All other bodv and cross product subroutines are patterned after these 
files. This apparently cumbersome organization had one significant advantage. Once 
the files were written and assembles. bodies could be toggled on and off in the 
computation simply bv adding a "°o" at the head of the appropnate subroutine call 
line. This allows the overall SLICE geometry to be reduced into smaller pieces rather 
easily and helped demonstrate the complex interactions that occur in the full 


multibodied hull. 


4.3 Results 

Only partial results were obtained for the numerical solution during the 
allocated period of study for this thesis. Complete validation of results was not 
obtained and the results contained in this section require further analysis before design 
decisions can be made. Sufficient progress was made to draw some preliminary 
conclusions. Model \{-2 was computed in the speed range of 3-9 fps (full scale 11-22 
knots) using the values shown in Table (4-1). The results presented in Figure (4-2) 
and Table (4-2) are based on using fifty subdivisions in the x integration, 12 divisions 
in the z integration and 10 divisions in the 6 integration. Testing was conducted 
varying the number of divisions and this was found to be a minimum acceptable limit 


with the integrity of the result maintained. 
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Fwd Strut Offset (v axis reference) + 0.979 ft 
, Aft Strut Set Back (x axis reference) 
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| Strut Submergence (z axis reference) 


|| Poa"Eereth 






| Pod Diameter (max) 


Pod Submergence (z axis reference) -0.59 ft 






Table 4-1: M-2 Gross Parameters Used For Calculation 


Figure 4-2 presents only the wavemaking portion of resistance as predicted by 
the numerical scheme. The negative overall result is attributed to the overall sign 
convention that is used in the calculation which has x positive in the direction of 
motion. An unresolved problem with this conclusion 1s that the body terms (1-1 


through 8-8) should have been negative on their own. This however was not the case. 
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Figure 4-2: Wave Making Resistance of SLICE calculated from Thin Ship Theory 
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Speed 3-3 4-4 5-5 6-6 7-7 8-8 | 
(fps) | (Ibf) 7 (ibf) | (ibf) | (ibf) | (bf) | (Ibf) | (lbf) | (bf) 


3 Oro E=3 2E-3 0.62E-3 0.62E-3 0.0038 0.0038 0.0038 0.0038 


oe eo 0.65E-3 | 0.0061 | 0.0061 | 0.0061 | 0.0061 | 
















Table 4+2a: Wave Making Resistance Component Terms 
BODY TERMS 


Table +2b: Wave Making Resistance Component Terms 
SAME CORNER CROSS TERMS 
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vSpecd(aia ets a 6-8 | 

| dps) | ab ab Abt) abf) | dbf) | dbp 
0.276 D Riher | 9 245 WP Ooo il) 0.421 |) 0.241) W046 | 

Ca 
0:356 ae -0.134 -0.528 0.350 0.383 -0.37] -0.003 


fA 


-0.388 ' bole 1.264 -0.006 -0.903 1.106 27 -0.004 


Table +2c: Wave Making Resistance Component Terms 
LONGITUDINAL CROSS TERMS 





Speed 1-4 1-8 5-4 5-8 2-3 2-7 6-3 6-7 


(fps) (Ibf) : (Ibf) (Ibf) (Ibf) (Ibf) (Ibf) (1bf) (bf) 


| > [ee 0.04E-3 | o119 | -0018 p00 002-3 | -o043 | 0018 | 00 
0.250 -0.574 0.001 -.LLE-3 -0.045 -0.57 0.001 


E-3 0.401 0.815 -0.003 -.30E-3 0.870 0.815 -0.004 | 















a: aes ia ania aia ae 
ae 0.75E3 | -1308 al EZEBP 0.006 | -31E3 | 0986 peel 0.005 


Table +2d: Wave Making Resistance Component Terms 
DIAGONAL CROSS TERMS 
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3-8 7-4 
(Ibf) | (Ibf) 


13 
=n 0.25E-3 | -0.492 0.003 0.131 0383 | 0.001 
ae 0.13E-3 | -0.484 0.003 | 0.28E-3 | -0.406 | -0.489 | 0.003 


-0.404 -0. 404 0.005 0.34E-3 -0.491 -0.944 0.004 


Table 4-2e: Wave Making Resistance Component Terms 
TRANSVERSE CROSS TERMS 





Understanding that the code used is only partially validated. the following 


observations are made: 


SrRUTS 

The contribution made bv the struts is small. This is true both of the body 
terms and the cross terms between any two struts. This is reasonable for the strut 
geometry used. There is symmetry about midships. the nose and tail sections are 
slender and the overall L’'D is 7.4. To verify that the low values obtained were not a 
code problem, an alternate geometry was selected that did not possess midships 
symmetry but retained the same L/D. A significant increase in resistance values was 
noted. These results follow the conclusions drawn by Chapman and Zhu Bing-quan on 


strut optimization cited in Chapter 2.5. 


ING So Miekivis 


As anticipated, the Same Corner Cross terms have the largest impact on wave 


78 


making resistance. This is followed in importance by the Longitudinal Cross Terms. 
The Diagonal and Transverse Cross Terms had more significant impact than was 
anticipated. Although smaller than the other terms. it would not be correct to eliminate 
the Transverse terms from optimizing iterations. Also clear from a review of Table 
(4-2) is the complexity of the bodv interactions. At a given speed. the transverse term 
for a specific body may have an additive effect while the longitudinal term for that 
same body has a subtracting term etc. Additionally. due to the lack of fore aft 
symmetry in the baseline SLICE layout. similar body groupings do not behave the 
same at a given speed. Note the difference between 3-8 and 7-4 in the Transverse 
Terms Table as an example of this. This decomposition of SLICE wavemaking 
resistance demonstrates that optimization work must be conducted with the entire 
hullform present. Attempting to sub-optimize a component body (such as a single strut 
pod combination) without observing the impact on the entire hullform mav lead to 


erroneous conclusions. 


a2 


Chapter 5 Conclusions 


5.1 Restatement of Thesis Objectives 


This thesis attempted to add to the developing body of SLICE knowledge by 
focusing on the resistance and powering aspects of the SLICE vessel in an 
academically mgorous manner. Three specific objectives were identified: 

1) Independent research into the SLICE hullform to give both the US Navy 
and LMSC a second hvdrodvnamic viewpoint of the 177 Lton vessel prior to the 
ATD. 

2) Development of a larger bodv of data to aid in optimizing the 177 Lton 
hull form. 

3) A broader understanding of the hydrodynamic features of the hull to 
facilitate deciding the memts of pursuing designs in other displacements more practical 


to the US Navy. 


In order to meet these objectives, several questions required answers: 

1) LMSC theorizes that the appropriate charactenstic length for SLICE 1s the 
length of a single pod. Is this valid? 

2) The specific interaction between bodies has not been defined by LMSC. 
What is the relative importance of form; strut pod alignment: and the effects that a 
strut/pod combination have on any of the other three strut/pod combinations that form 
the underwater body? Understanding these effects is necessary for design optimization. 

3) The LMSC proposed design places the forward struts inboard of the aft 
struts. Is the outward lateral force (y axis lift) created by an effective angle of attack 
on the aft struts significant enough to create structural concerns? 

4) The LMSC model tests were all conducted with models fixed in pitch and 
heave. Is this an accurate enough representation of realitv to make predictions? 


5) The LMSC tests maintained speed increments for an extremely short 
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duration (usually less than 10 seconds) and frequently truncated test runs at a model 
scale equivalent of 25 knots full scale. Was this sufficient to characterize a ship with 
an intended speed of 30+ knots? 

6) Is the apparent focus on skin friction reduction. as reflected in the body 
forms selected. reasonable? 


7) Does SLICE offer the stated hydrodynamic advantage over a SWATH ship. 


5.2 Conclusions 


CHARACTERISTIC LENGTH 

The appropnate charactenstic length for SLICE is the strut length. not the pod 
length. This is illustrated in the resistance curves for M-1.M-2 and M-3. These models 
had significantly different underwater bodies but very similar strut lengths. On a 
resistance vs speed plot. the position of the wave making hump for all three hullforms 
is nearly identical. In order to align these with .Y= 4.5. strut length must be chosen as 


the characteristic length used in the definition F. 


THE IMPORTANCE OF BODY FORM AND BODY INTERACTION 

Model \{-2 was developed with the expectation that finer lines would help 
improve the hump speed performance of the baseline SLICE hullform. The offsets 
selected for the pods increased the constructive L/D ot the pods from 4.1 in the 
baseline to 5.0 in M-2. This was the maximum improvement that could be made to 
L/D and still maintain the gross characteristics of the Baseline SLICE. This change in 
L/D was not sufficient to impact the hump speed performance. If the SLICE hullform 
is maintained. improvements in hump speed resistance must come from body 
interaction. Underwater photography showed that propeller submergence 1s better 
through the hump speed than previously thought. 

The finer lines of M-2 did make a significant improvement in the resistance 
measured in the high speed range. Underwater photographv revealed that the bluff 
lines of M-1 causes the shedding of a significant vortex at the trailing edge of the 
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struts starting at model scale equivalents of 25 knots. These vortices were incorrectly 
interpreted as spray phenomenon in previous testing. In the full scale ATD ship, the 
flow modification caused by the propellers may cause these vortices to be ingested 


into the propellers with a significant loss of thrust. 


MODEL TEST CONDITIONS 

Attempting to test the SLICE hullform in a condition other than fixed in pitch 
and heave introduces a significant level of complexity to the model design. The added 
complexity does not necessarily make the results more useful. The difficulty lies in 
developing a system that connects the model to the tow carnage through the line of 
force. With four underwater bodies all towed from a central point. this is extremely 
complex to manufacture. Anv connection point somewhere other than in the line of 
force 1s easier to manufacture but creates variable applied moment. This effect cannot 
simply be subtracted out since it is impossible to separate it from the effects of real 
hvdrodynamic forces. The remaining option then 1s to lock the model in pitch and 
heave. For the purposes of conducting comparative testing, this simplification 1s 
justified and was a reasonable simplification made by LMSC. 

The tests conducted as part of this thesis showed that there was a significant 
transient time associated with acceleration from rest. In the higher speed ranges, forces 
required 4-6 seconds to achieve steady state after the steady state speed was achieved. 
The LMSC tests incorporated more speeds in a single run with fewer accelerations 
from rest. LMSC allowed for significantly less data collection time at each speed 
increment (less than 10 seconds compared with greater than 20 seconds). This may not 
have been sufficient time to achieve steady state and may contribute to the differences 
in LMSC's test results and those conducted as part of this thesis. 

The numencal analvsis conducted in Chapter + highlighted the extreme 
complexity of the flow around SLICE. Relationships presumed negligible such as the 
transverse and diagonal body interactions were shown to contnbute noticibly to the 
overall resistance value. This has an implication in the LMSC test results. When the 


August control surface tests were conducted. only one side of the model (starboard) 
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was adjusted for various angles of attack. Based on the results of the wave making 
decomposition however, drawing conclusions about whole body behavior based on 


partial body measurements may lead to incorrect results. 


HULLFORM OPTIMIZATION 

The data decomposition presented at the end of Chapter 3 demonstrates that it 
is incorrect to optimize the hull for high speed operation by focusing on frictional 
resistance. Powering through the wave making hump does not imply that wave making 
resistance becomes negligible. It is true that the coefficient of wave making resistance 
C,, becomes a small number but resistance continues to increase as velocity squared. 
The design perspective should be that Froude scalable effects will account for at least 
half the total resistance and the hull should be optimized accordingly. 


SLICE APPLICABILITY TO WHE Us NAV 

Figure (5-1), obtained from a set of presentation graphics used by Dr. R. 
Compton at the U.S. Naval Academy, depicts a predecessor to SLICE. The literature 
survey conducted as part of this thesis did not uncover any design work accomplished 
on this hull. 

The comparison of either M-1 or M-2 with M-3 indicate that if there is an 
advantage to the SLICE hullform over conventional SWATH hullforms, it must be 
found somewhere other than in powering performance. M-3 demonstrates that it 1s 
possible to create a SWATH to the design constraints of SLICE and still meet or 
exceeded SLICE resistance performance through the entire speed range. 

SLICE offers a possible advantage in the ratio of Deck Area to Displacement. 
The development of M-3 showed that the 177 Lton SWATH could have up to 600 ff 
less deck area than the equivalent SLICE. The extent to which the crossbox is 
cantilevered fore and aft of the struts becomes the controlling factor. For low density 
pavloads. such as ferry traffic, this is significant. For higher density payloads typical 
of military craft. the advantage is diminished. 
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Figure 5-1: A SLICE Predecessor 
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The following recommendations are made for future study of SLICE 
technology: 

1) The question of lateral lift on the aft struts created by the flow direction 
remains unanswered. This question should be addressed before the ATD vessel is 
placed in service to ensure a structural issue, either from an integrity or fatigue 
standpoint, 1s not overlooked. 

2) This thesis presented a reasoned approach to scaling model data to full 
scale. The ATD will be the validating test of assumptions and decisions made in 
presenting the case for using traditional methods of data scaling. If conducted 
properly, the ATD will be invaluable in demonstrating the robustness of Froude 
hypothesis. During the ATD, significant effort should be invested in photographically 
documenting flow behavior both above and below the surface. Without this 
observation of physical phenomena. specific conclusions about scaling issues will not 
be possible. 

3) The numencal approach used in this thesis requires further development 
before tuning of the SLICE hullform can be fully demonstrated. Thin ship theory 
appears to be reasonable for comparative analysis purposes and this approach should 
be further developed. After further validation, SLICE component positioning can be 
optimized and then a model of the final outcome should be built and tested in a 
similar scale to M-2 so that results can be compared. This could be the subject of a 
reasonable follow on thesis. 

4) The SLICE ATD has the potential to answer many questions that have been 
raised regarding the classical naval architecture approach to resistance predictions. 
With careful planning, the ATD will be much more than a demonstration of a single 
ship's performance. It has the potential to contribute to the Naval Architecture 


community's understanding of the robustness of current testing and scaling techniques. 
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Appendix A 


Model Testing Lessons Learned 
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The Science and Art of Model Testing 


The purpose of this appendix is to outline the thought process that was used in 
designing the experiments and the models that were used in this thesis. The tone of 
this section will be conversational and it represents a compilation of lessons learned as 
the research for this paper progressed. It is hoped that this appendix mav be of use to 


other students who may choose to undertake a similar type of research effort. 


When the work for this paper began, it was decided that in support of a Naval 
Construction program, there would be a wealth of experience to be gained from 
designing and building a ship (albeit a small one) that actually had to work. With this 
decided. I embarked on a challenging. frustrating. rewarding, extremely time 
consuming, detail intensive process. 

The title to this appendix is not meant to be tnte. There are clearly 
scientific/engineering principles that may not be violated if reasonable test results are 
expected. Geometric similitude must be maintained. tolerances must be decided on, 
load and matenal strength estimates must be reasonably predicted, scaling must be 
calculated properly and data reduction must be interpreted accurately. Bevond this 
however, I discovered an entire world of additional issues regarding the planning of 
the expenment. actual construction of the models and successful time and resource 
management during the cnitical and alwavs too short time allocated in the test facility. 
You cannot gain an appreciation for this facet of experimentation bv taking a course 


or reading a book. 


RECOMMENDATION #1: Prior to doing anvthing on your own expenment, 
observe similar projects for several davs. In the case of this research, I had the 
opportunity to observe the SLICE testing conducted at NSWC in August and 
December. For eight days that I was able to watch other people devise solutions and 
cope with frustrations without impacting my own work. Many ideas for my own 
project were born during these eight davs and most importantly, I met the experts in 
the field of model testing. Later, these associations Were invaluable to me. 
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If this penod of observation does not scare you awav from proceeding with 
vour own experiment. the next question vou must honestlv answer is: How handv are 
vou and what resources do you have at your disposal? When I began building the 
models. I considered myself reasonably proficient with tools but I had never 


undertaken a project of this magnitude before. 


RECOMMENDATION #2: As a means of determining your handiness, I 
suggest vou consider the type tolerances you will have to build to uf the model is to be 
considered a reasonable geosim. For argument. Jets sav that in full scale, the greatest 
building error that will be tolerated is 6 inches (pretty big error!). If the scale that you 
are building to is 16:1. inaccuracies in vour model of greater than 38 inch exceed the 
building tolerance and vou have introduced a level of uncertainty as to the actual 
geometnc similarity of vour model to the real ship. If vour personality is such that vou 
fail to recognize the difference between 3/8 and 1 2 inch when measuming things, then 
model building mav not be the thing for vou. Scaling errors to full scale to determine 
reasonableness was a process that slowlv evolved during mv own project 


At this stage. some decisions need to be made regarding the specific test 
objectives. To a large measure. this will determine the complexity of the model to be 
built. The speeds to be modeled must be matched with tow tank capabilities. For this 
project. full scale speeds of 40 knots were desired. Based on the maximum camage 
speed at the {IT tow tank. this would have required building the models in a 100:1 
scale. a size so small that no useful data would be collected. The alternatives were to 
either slow the maximum speed or find another facilitv. The United States Naval 
Academy was literally the closest facility with the necessary capability. In addition to 
these concerns. there were also the questions regarding what the model should be 
capable of doing. Early on in the project it was decided that the models would have 
no contro] surfaces. This removed a significant level of complexity from the model 
design. It was also decided that the models should be capable of being modified for 
self propulsion at some point in the future. Additionally, since the general test scheme 
called for a test free in heave. the models had to be built with sufficient reserve 
buovancy to support themselves, the mounting frame. the towing carnage heave post 


and force blocks and an estimated reserve for propulsion equipment. This simple 
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exercise placed manv constraints on the model: 

1) The models would have to be particularly rugged to enhance the chance of 
survival during the 400+ mile tnp to the test facility in the back of a van. I eliminated 
the conventional foam construction in deference to this requirement. 

2) By understanding what the models had to do. I could make a weight budget 
to build to. It turned out that the heave post and force block weighed 25 pounds 
collectivelv. Properly scaled. I could estimate the displacement of the models (between 


75, 83 and 115 Ibf for models \f-1, M-2 and M-3 respectively). I then budgeted 10% 





for reserve buovancy and allocated 15 Ibf for frame weight. The remainder could be 
used in materials to actually build the models. This weight budget further drove the 
building plan. The material selected had to be strong, extremely light and simple 


enough for an amateur to work with. 


RECOMMENDATION #3: Design vour models only for the capability you 
plan to test. Build only what vou need. Clearly establish exactly what the models wil 
be used for before planning their construction. Eliminate multiple functions where 
possible since multiple functions will most likelv increase the model's level of 
complexity. In the case of this project, significant effort was expended to meet an 
established weight budget. The weight budget was based on the assumption that the 
model would be tested in a free to heave condition but this test was never performed. 
For the fixed 1n heave tests, the models would have been acceptable even if they were 
negatively buovant. 


At this stage there is almost enough preliminary groundwork set to begin 
building. With anv luck vou will be able to enlist the help of others in the actual 
building process. The models built for this project represent the collective work of 9 
individuals with a total of roughly 700 manhours invested between December and 
March. Admittedly a portion of this time was expended inefficiently, devising 
solutions to emergent problems. obtaining supplies etc. but these are the facts of model 


building and this sort of time expenditure should be planned for. 


RECOMMENDATION #4 If vou are not literate in a CAD svstem, develop 
that proficiency well in advance of the model building. As early as possible, make full 
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scale working drawings. This is critical if other people are helping you. It is the only 
wav to know that vou helpers understand what to build. If the drawings are done soon 
enough in advance. they mav help point out errors in the design. I did not develop a 
good set of working drawings early enough in the process and it hindered progress on 
several occasions. Once the drawing process was brought under control, I found it 
handy to cut up the rather disposable drawings for templates. 


Oniginally P\'C was the material of choice for the model construction. It is 
easv to shape and bond and 1s a relatively inexpensive alternative. Additionally, water 
absorption is not a concern. There are two drawbacks though. PVC is relatively dense. 
At 81 Ib/in’, it can be made to float if rather thin walled sections are used. This is 
useful for long stretches of parallel midbodv but becomes a significant drawback for 
the shaped pieces which must be turned on a lathe. Hollowing these sections out to an 
appropriate wall thickness was not practical. Additionally. the time to tum a single 
nose cone (inside and outside) was fund to be about 10 hours after considerable 
practice. A total of 16 nose and tail cones were required for two models. 160 hours 
was considered excessive to finish this single piece of the project. The second 
drawback of PVC 1s also one of its advantages. PVC is relatively easy to cut because 
it is soft. Although there was little danger of it deforming in a gross sense, it 1s very 
easv to gouge and scratch. Once damaged in this fashion, it is difficult to repair 
because although it is easv to sand and it bonds well to itself with the proper adhesive, 
other bonding and filling agents do not adhere well. In short the probability of 
achieving a satisfactory repair once damaged is small. In light of the tnp the models 
would take to the test facility and the potential for damage during this tnp, this was 


considered a major drawback. 


RECOMMENDATION #5 Treat your model building project like an 
actual shipbuilding project in the preliminary phase. Account for weight and buoyancy 
to the greatest accuracy possible. I developed a weight budget for Structure (the 
frame), Propulsion (the tow pad), the Pavload (the heave post), a building margin of 
10° and a design reserve buoyancy of 10% (allocated because there is some 
intention to eventually convert my models to self propulsion so there may be some 
growth). By calculating the displacement. J then knew how much each component 
could weigh and I did weigh them frequently through the building process. As work 
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proceeded, building margin was allowed to be utilized at a rate commensurate with the 
stage of building. I incorporated this method early enough in the project to abandon 
PVC as a building matenal before too much work was done, but not soon enough to 
avoid buving a quantity of matenal that was not useful. 


RECOMMENDATION #6 Make frequent tests of methods vou think you 
will use, well in advance of when vou need them in the building project. If you plan 
to use a particular machine, check it out to ensure all the pieces are there and it 
functions properly. If you plan to use a particular adhesive to bond two different 
materials. try it out early on with two scrap pieces. even if the adhesive claims to be 
formulated to do the job. I had significant difficulty with this on several occasions. 
Testing gave me the time to either call the manufacturer for assistance or to find a 
different product. The two times I did not perform tests on intended procedures 
resulted in near disaster. The one time I did not perform a bonding test was with PVC 
and PVC cement. items I had worked with in the past and was assured by the PVC 
supplier would work. For a reason I never determined, the cement did not bond. This 
occurred late on a Saturdav immediately proceeding the Sundav I was to depart for 
Annapolis. The search for a new bonding agent was frustrating since it was not 
supposed to happen. difficult because of the hour. and nearlv cost me a significant 
piece of the test program (the high speed SWATH model did not require significant 
reserve buovancy and so was built from PVC to utilize the matenal purchased.). I 
intended to paint the models bnght vellow with marine paint. Yellow is the color of 
choice for model testing because it enhances contrast and is useful for photography. 
The marine paint was difficult to apply. did not cover well and was slow to dry. This 
provided for extreme frustration at the alreadv franuc end of the building phase. In the 
end analvsis. anv sprav enamel with a 10-15 minute dry time would have performed 
equally well in the tow tank tests. and I would have known that if I tested the paint 
systems in February rather than waiting until the end of the project. 


For the SLICE hullforms. molded fiberglass was tound to be the best 
construction method. With other materials such as shaped foam. PVC rod turned on a 
lathe or stacked wood, construction is a 100% hands on process. Only one piece can 
be manufactured at a time. Each SLICE model has the following components: 4 struts, 
4 nose cones. 4 tail cones and 4 parallel midbody sections. There were two SLICE 
models, therefore requiring the manufacture of 32 components. Hand crafting each 
component would not only have been prohibitively time consuming but also would 
have introduced slight variations between components that should have been identical 
(the four nose cones for example). Molded fiberglass construction offered the 


following advantages: 
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1. A mold insures that the four components have the same shape. 

2. Once the molds are made. component fabrication time is pnmarily 
epoxy cunng time. This frees up time for the builder to do something else. 

3. By laving in the fiberglass cloth in the correct fashion, excellent 
strength is obtained for very little weight. 

4. Using the proper additives in the epoxy provides an outside surface 


that requires no finish work once the piece is removed from the mold. 


The West Epoxy Svstem was the matenal of choice for this project. It is 
relatively expensive but is readily available at Marine Supply stores, very easy to use, 
and 1s associated with an excellent consumer product hotline that is responsive and 
helpful. Most of the products referred to in the building process are West System 
names. The West svstem can be used with either a fast setting hardener or a regular 
hardener. The fast hardener cuts the working time and also the curing time of the 
mixture in half. Additives which are called for in some of the steps also reduce 
working time. Times cited below are based on using the fast hardener. The West 
System sells a set of proportioning pumps which helps to ensure that the correct mix 
of resin and hardener is achieved. The pumps are a worthwhile investment. 

Note: The epoxy will not adhere well to plastic so working on a plastic drop 
cloth will help cleanup. Mixing the epoxy in plastic pails can be reused. Just let the 


epoxy harden and the tap the pail. 


The following outlines the steps used constructing the pods: 


1. A CAD drawing was made of the pod. Five or six tangent lines were drawn 
along the curvature. Angles with the vertical and insets from the leading edge were 
calculated. This is shown 1n figure A-1. 

2. A male plug was cut out of PVC rod on a lathe. This was done by marking 
the insets. and then cutting at the appropriate angle from the inset mark to the leading 


edge. The transitions between angles was fared with a file and the entire plug polished 


a2 


with consecutive applications of 200, 320 and 400 gnt wet sandpaper followed by an 
abrasive cleanser and finally brass polish. All these steps were performed on the lathe 
to ensure roundness. 

3. A female mold was then made from the male plug. This was done with the 
West Svstem. First the Plug is coated with a laver of wax. This 1s necessary to ensure 
the epoxy does not adhere to the plug. The first epoxy coat consists of the appropriate 
mix resin and hardener with sufficient microfibers added to give the mixture the 
consistency of paste. Since the additives decrease the working time. I suggest you 
pump the resin into the mixing pail, add the additive and mix it thoroughly and then 
as a last step, add the correct amount of hardener. This method will save a few 
minutes of working life. This mixture is painted onto the plug and serves as very 
smooth gel-coat laver. Brushes can be conserved by washing them with acetone after 
use but I recommend alwavs doing this first layer with a new brush. This application 
is allowed to begin to harden but not cure (approx 2-3 hours). Even with steep angles. 
this mixture should be thick enough so there is very little running of the coat while it 
sets. A second laver is added which 1s also a paste like consistency. This time, the 
additive is a 50°o-50°% mix of microfibers and colloidal silica, which will be less 
smooth but has greater strength. This also is set aside to begin to harden but not cure 
(2-3 hours). Finally, the real strength layer is added. For the molds. this consisted of 
squares of fiber matting soaked in mixed epoxy (no additives in the mixture this time). 
The matting has a random fiber lav and tends to be fairly (0.125" thick). Do not let 
the matting get too wet or else it will tend to slide when placed on the plug. The most 
effective way to do this step is with your hands. Latex gloves are an absolute 
requirement. This step is just like working with paper mache. When the matting has 
begun to harden (2-3 hours) a laver of 12 ounce woven fiberglass cloth is added as a 
final step. Pre-cutting the cloth into smaller patches will help since vour gloved hands 
will be covered in epoxy. Since there are no additives in the epoxv in the last two 
steps, use the epoxy sparingly, just enough to barely wet the cloth. Allow the entire 
mold to cure for the appropnate time ( 8 hours with fast hardener). The mold may 


have to be beaten pretty hard to pull it free. In the case of my shapes. I found that 
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removing the plug actually drew a significant vacuum because of the closeness of the 
fit. The mold is strong enough to be beaten prettv hard with a rubber mallet. A small 
hole dniled into the mold and later patched may also help. 

4. Manufacturing pieces is now essentially a repeat of step 3 with a few 
exceptions. I chose to leave out the fiberglass matting in the final pieces. There was 
sufficient strength with the two initial coats and the 12 ounce cloth and I wanted to 
reduce weight. Take Note: The epoxy cures using an exothermic process. In the 
steps used up till now, this was insignificant. Since the female mold is a cup, as you 
manufacture the component pieces. anv excess resin will pool in the bottom of the 
mold. On two occasions I cracked my molds because of the heat released by the 
excess pool. I recommend that vou set the mold into a pot of cold water at least 
during the final step (glass cloth and epoxv). This seems to be adequate to control the 
temperature. Jake Note: You need to provide a method of pulling the piece out of 
the mold. I did this by setting two 12” pieces of nylon string into the wet epoxy of the 
final coat and covering the end with a square of fiberglass cloth. This gave a secure 
method of pulling the piece out. 

5. After the components are pulled from the mold. thev should be washed 
with acetone to remove any residual wax. If this is not done later epoxving during 
assembly may not bond well. 

6. When all the components are manufactured, thev can easily be assembled by 
epoxving them together. I cut rings out of balsa wood that were fitted into the joints to 
increase the gluing surface. The timing of this assembly is dependent on how the 
struts get attached. For the SLICE models, the struts are bolted to the pods. In order to 
facilitate this, the tails were joined to the parallel midbody sections. The struts were 


then attached to the sub assembly and finally the nose cones sealed the assembly. 


The following process was used to build the struts. 


1. A full scale CAD drawing of the strut cross section was made. This was cut 


out and glued to a piece of thin plywood. The cross section shape was cut out of the 
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wood so that a “slip over template was created. 

2. Blocks of balsa wood were cut to rough dimensions of the strut and glued 
together to create the correct strut length. 

3. Once dry, the balsa wood was sanded to shape so that the template made a 
loose fit over the balsa wood. 

4. The balsa wood was sealed using a coating of epoxy painted onto it. This 
was left to cure. 

5. The strut was sanded to remove large defects. (The balsa wood will cause a 


fair amount of bubbling in the epoxy). It is then covered in 12 ounce fiberglass cloth. 


oe) 











Figure A-1: Example of Nose Cone Template 
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Appendix B 


Model M-1 Baseline 
Table of Offsets 
and 


Design Worksheet 
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Fowara Pod in Inches 


(a 


rl | pes 
tmp Gilt: 
Ds aes 
oc} 


Mode! M-1 Table of Offsets 


Aft Pod in Inches 


aX 


api ap, 
0 














n 
O 
we) 
N) 
oO 
JN 
NS) 
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toftro 
my 
ea) 


tU] top to 
42] OV 
a on 
—] vs 









0.695 
1.39] 
2.086 
27S! 
3.456 
9.06 
oie 
10.43 | 
11.125 
hee 
12.516| 
iB 
13 906| 
14.602| 
15.297 
15.992| 
16.687 | 


een Bann ee See 


Struts in Inches 








These are the full scale SLICE offsets taken from the December 1994 LMSC 
variant used to generate Model M-1 offsets. 





Pod X Coord (ft) 





——_—_— ee — 
Pod Y Coord (ft) | Strut X Coord (ft) | Strut Y Coord (ft) | 





Model M-1 Baseline LMSC SLICE 


This is based on the drawings provided by LMSC for the December baseline ship. 


1a. Units Defined 
ft=1L Ih= 1M sec=17 Base Units 
nm =6076-f Nautical Mile 


Kt 24 Knot 





nm 
lur 


lton =2240 1b Long Ton 


1b. Foward Pod Offsets in Full Scale 


\ tp. = READ( mI pix) \ tp. = READ (ini piv) 


Aft Pod Offsets in Full Scale 


“ = READ( mI pan ) 


ap. “ ap. =READ(mI pay) 


Strut Offsets in Full Scale 


Soe se (11 ESN) Vv. =READ(mI sv) 
j f 


Model Scale Factor 


: Soi ; ; 
A = = feet by Boe te: 
SBC ah 
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Model Offsets 


“ fo tp) ion L2 


101 


7 ae 
ioe Ip / he NA 


. at eae 
pi? ap/ Baal) /. [2 


1 


Ya 


Foward Pod 


Aft Pod 


Struts 


1c. Model Scale Weight Balance 


Fixed Variables 














: 
SUG a 
ae aoe. ae Density cf Sait Water at 56° F An assumed average 
re cperating temperature. 
ple 
N44 5 it c 
Veter aod 10 ninematic ..scosity of Salt Vater at SE° F 
BiG 
ec : 
#4 D507 a0 ~ensitvy ct =-esn Water at 66° © = assumed average tow 
i" "24K temcerature Use tnis if crea:cion is for a moaei test. 
,e 
pemcese TLL 
Ae QSOS 10 ~.nematic + scosity of Fresn Vater at 68° F 
SCC 
. = a . 
g=32 i 74sttsec ~ Acceleration due to gravity 
a V =O at. i kt.. 35-kt Modify this to reflect the velocity range of interest. 
Indicator Settings 
Water 
SW "1° inaicates the sn:p is operating in Sait Water 
"O° inaicates the snip is operating in Fresn VVater 
= ap & i Z 
op =a SW] He coy he 
v =iSWeHl.y 


Soe 
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Hull Displacement Definition 


_ 170-Iton This defines the the requireded (q) displacement of the vessel 





d R 
aS note: 170 Iton its usec here since that ts the displacement of the “filletiess” hull. 


Vg ee Wa ear The requireded (q) displaced volume 


pe 


Bouyant Force 


4.084-1b 


~) 


Ph FV gee ey 
Weight Account 


FRAINIE 


e's ss - IS Ib 
Force Block - Heave Post 


W b] = 2>1b 


Other Mode! Supplies 


\W oth = 2-)h 


Desired Reserve Bouyancy 


Allowed Veight of a single strut pod 


ee a Wag 8 
Sp 4 


oes nnd 


\W bey a 1Gelb 
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Appendix C 


Mode! M-2 Modified Baseline 
Table of Offsets 
and 


Design Worksheet 
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Model M-2 Table of Offsets 


Pods (both fwd and aft have the same offsets) struts 
in inches in inches 


| yscale(x1) 


: : =i 
vm(v) in Nm(N). 1n vmev). mn 
yy} (N) - ee \scale(x1)-- | 





lyme 


to 
{.) 
04 
‘7 


See 

' 

5 | 

8 } 

nnd 
i 


EE 


Sli 
re 
cm 

—] bo 


to] to 
tJ] t.2 
tystt > 


ror] t 
a ot 
Lofts 


wlol— 





10 569} 
il 125] 
ae 


2G 1a] 
112 238] 26.7 (0 GNI 
te fd ae 0 = 
e55) LI pay ee ) 
13.906} eas 





l 


1.909} 
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Viodel M-2 The Better Slice 


This vanant uses the same gross dimmensions as the baseline SLICE but varies the 
body shapes to follow more traditional forms. The struts use parabolic leading and 
trailing edges and the pods are based on series 58 shapes. 


1a. Units Defined 


ft=IL Ib= 1M sec=1T Base Units 
nm =6076-f Nautical Mile 


= 2 


kt Knot 


lton =2240-lb Long Ton 


1b. Fixed Variables 





Sec 

San = | 9912-Ib-—— Density of Salt Water at 56° F. An assumed average 

a‘ operating temperature. 
5 it 

epi SE ee LO + — Kinematic Viscosity of Salt Water at 56° F. 

sec 
Be 
P sec ; 

Brae taco Ib: Density of Fresh Water at 68° F. An assumed average tow 
ft tank temperature. Use this if prediction ts for a model test 
aie 

Oe = |.0804 1Q ~-—— Kinematic Viscosity of Fresh Water at 68° F. 

sec 
-2 
g = 32.174*ft-sec Acceleration due to gravity 


cia V -=O0-kt, 1 -kt.. 35-kt Modify this to reflect the velocity range of interest. 
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1c. Indicator Settings 


Water 
$$$ SW =! "1° indicates the ship is operating in Salt Water 
°O” indicates the snip is operating in Fresh Water 
p= if(SW>1 p SEHD) 
v =if/SWEl.v av ay] 
Body Segmentation 


$$$ seg =1-ft 


2a. Hull Component and Required Displacement Definition 


Xa = eiton This defines the the requireded (q) displacement of the vessel 


V gird i ve =6189+ft° The requireded (4) displaced volume 
Strut Height 

H =6-f Ammount of Strut in the Water Defined per the Baseline 
Strut Beam 


vw" DD, 23.25-tt Maximum Strut Beam Required Defind per the Baseline 
Considered a required parameter. 


Overall Strut Length 
Pee eect Overall Strut Length 


Length to Beam Ratio 


E 
ED . LD! =7.385 Length to Maximum Diameter Ratio 
l 


Length Breakdown 


one L omb| = 7-ft Length of Parallel Midbody 
one Ea = (>: (L 1- Lomb! ' Fe ia 8.5°ft Strut Nose Length (2.4 x D is optimum) 
eres ee Lomb] -L, Ly, =8.5-ft Strut Tail Length (3.6 x D is optimum) 


Nose and Tail Shape Exponent 


as iV) eae iia Ny) 22.25 


Calculation of Nose and Midbody Shape 


eee et all amu defines the x values of the subdivisions 
L ra ies 
Ave ad & ee a | defines the y values of the nose based 
ore i > RS | on a parabola 
nl 
D | 
Von sl). == —— defines the y values of the midsection 
ESTER aa eS Bl Be pes) pmbliGn |) 
Calculation of Tail Shape 
x1 =O-ft,-seg..- (Ly - L mbt ~ Lay) defines the x values of the subdivisions 
( ia Mey | 
ala ae x1] ~ Lait lL pmb! | defines the y values based 
. 2 cy ho on a parabola 


bodyl(xl) =aff [x1] Ly) +L pmpy AL(S!).stI(x1)) Function that defines the 


pmb 
body geometry 


LOS 


Strut Section 





Wetted Surface and Volume Calculation 





Perimeter 
Nose 
| 2 2' 
D a Dy = ome 
N05: Py 161g Pe oot tg) = ar 
8: LF ‘ D ] 
Np =17.55-f 
Tail OD? 4b yt DY = 162, 
lee = 05: Dae 16L 4) So aN ee 
P \ 8-L D 
tl ] 
qT, = 17.55+tt 
Midsection 
Sp = 2-L omb| 
Sp = |4-ft 
Tota! Perimeter 
P =NorSprl, 
P =49.1-ft 
Wetted Surface of the Strut 
WS, =PH WS, =294.602-11" for one strut 


1Qy 


Areas 


nose area 


N. = 


9 oo ey OT 


Ww | rho 


tail area 


2 
Ty Faby Dy, 


led 


mid section area 


Sq =D, Lomb 


Total Area 


Area aire lq 


Volume per Strut 


Voll = Area A 
Total Strut Volume 
VO =4 vol 


Pod Volume Required 





V pods a d- VOL 
V : Vv pods 
pod ~ yi 


N , =18.417-f 
T = 18.417-f° 
Syn oie 


Area = 59,583°ft 


vol = 357.5+ft° 
VOL = 1430-ft° 
V pods = 4759+ ft" for four pods 
3 
V pod = SO it per pod 


Pod Defintion based on the Mode! 4155 of Series 58 


i l= 43 the scale factor used to get from the data on page 49 of senes 58 (Lpod= 9 
ft) to the size needed per baseline Lpod= 38.7 ft) is 4.3 


Ly =9-fta 


a= Olle SS 


mS READ (hbetx) 


xfs( x) =xL 


WS 6g =39.75-f° 


D, =L8ti, Dy, =7.74f 


y, =READ(bety) — non dimmensional 


yfs(y) =y-D p full scale dimmensions 
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a ee ch ee ae 
WS ps = VS seh | Ws go CS rent for one pod 


Vols. =14 9. As printed in the Series 58 data 


Voly, =Volsgi, ‘Voll, = 1183859: 
VOLUN Is aerev Ol weeracc! 


VOLUME =6165.437> ft V ae 6188.727°f Compares calculated volume 
to needed volume. 


4 =VOLUME p g 4 = 176.334 -lton 


Now Scaie to [oae! Size 





5 502540 ; | 
ia = ee a SG GR 
7 » p ‘7 A 
yy 2 Oi es 26.27 20 
AMUN) = NISCN)- A 5 vm(y) =stsCy)-4 5 
Better Slice Strut Offsets 4 = 006 
a ne eae ; 
\seaie 1) = -——-—— vscale(N]). = boda (5) | 2a 
ti i. 


11] 





Weight / Volume Balance 


Density of PVC 


[) : Orbs 75.111 ee SOS Sats 
eae ; 
aie Se Ee V pees ott 
4 
eer _ 
Prive = a P ove =X4 708°; 
r 
ee a Ibert 


Displaced Volume 


~ 


Vy rer e= NOIUNIS eae 


mode 


4 
eae 


! = 
\ model | 


Bouyant Force 


\? ie 


" v are ae 
a Made, LEN Ne 


Weight Account 


FRAME 


we ~ 0 


Force Block + Heave Fost 


diameter and length of rod 
weighed 


volume of 1 foot of rod 


ott” density of the PVC 
density of FW at 68°F. 
ye) 


Parallel IMidbody X 4 pods 


Hamp =Go In 
Da psycices > ee eee : 2 
S ee ee 5.5625 - 0.35 Z x 
“ 3veonin aa © ae rea 1 z re pmb 
, = ans 
\ pvepmb O.0G5°n 
: = me ) ee 2-15 
” pmo M pyepmb ? pve MW pmb 5-4B2-ik 
Pod Nose and Tailx 4 poas 
i: 7 eo See 
“sount > neds #2 74 L smb : ae volume of the nose 
2 and tail 
ee ie 
, SMUeak 
ae a 5 ae Aaa 
vs i eee ai ee 3976] in 
\ Se oO ne to 5s ca 
aes 2 OF oan! V ARO 133 in 
4 te : ~ 2 = 8 - a 5 2 = 7S my ead 8 
eoreae ad ees oie, Wane V anenes 28 2 ain 
2 3 
Ae sa SS 2 Sinn 2 =9R17 ein 
Snes oot 2S in ania V ahared 9X17 ein 
Se See ore 
re OP tsa. ais ee shoiesae Soon 
sing - y CN ES 2 ae as) = c7: we 3 
> see) ey Tl et Fe Mare 63 6]7 en 
= P) . - ei co = 53 * 3 
V ‘Rape ot C220. See V fan 75>. 578em 
’ . 2 eee ale eo ae 
\ ‘hopes. ee Seine 2 in 2 V inotee 282 74-in 
Ms - Xt: 4 2 = 17 S646. 3 
teeta aa Orie 2 ane vi (pee 1? 566 *1n 
As - = 5 - } 2 —— s 7 los 3 
Seams OS i a V eee eS? lena 


Vibore  Y nborel ~ ¥ nbore2 ~ Y nbore3 7 V nbored ~ © nbores 
* V shore2 ~ ¥ thore3 ~ Y thored 7 V thores 

Vater 280957 =n 

V pvepod = ¥ nodnt~ Y bore 

V pvcpod =O 479 

a podnt ~ : svepod P nye W podnt — Ne Ee 

4 Struts 
— TAAPea s 4 On V ia, sore 
eee at OG are er ream 2S 
ean Heian OME ne V shored = 0530 
" pvesir 7 Sa sbore~ “ <bore2 Y pvesir O Tolan 
W struts 7 " pvestr P pve 
Ne Ore) alle 
Other f.lodel Supplies 
Selma 
Checking for a Balance 
F irq 7 10lb Desired Reserve Bouyancy 
Ae } se Bs 

emoceh = tah pmo W podnt ~ W struts 7 Y oth 
ee efi Ca Sees Actual Model Weight 


Fira 7h h- Winode! Actual Reserve Bouyancy (Should 


F tyeq 2719 OOS lb 


ee ee bust pe Ulon greater, 


Claes rd 
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a tbore | 


= Insc) 


Areas to Help the Balance 


reserve = areas with excess weignt included that can be ttmmed down by the 
indicated ammount 


re =2lb frame ry, =9 1b force block +7, =1-lb other model items 


reserve -T eae Y eile r oth 


fixed = Wee Wyy~ W 5, - reserve Weights that are unchangable 
excess = reserve Douyancy that is desired but not needed 


excess =') Jb 


_ 


Po Cie CN Oes. 
a} I} ‘rd SENS 


Calculate the upper limits for strut‘nose and tail cone weights 


Wmodelmay =! b> F br 
a _ modelmay > HNed The Max Allowed weight for a 
spa 4 single strut/pod combination 
\W spa = 7.63 1+lb 
Wey Woe 
a se ec Seles Calculated Weight of a single 
4 4 4 Strut/pod combination 
Summation of Allowed Weights 
We eosin me i eclih Frame Weight 
Wo unpeee Force Block and Heave Post 
WN payed ieee Miscalleanous Weight Allowed 
W spa 7 7 631+Ib The Max Allowed weight for a 


single strut/pod combination 


Miode!l Dimmensions 
| 
hm = 
16.68 





a - oy ] = 20-f 
Lstrut _ Ly feces Lstrut .., 1. 439-ft 


WSsuut = WS oho” WSstrut ., = 1.059+tt? 


EpoG eee pod ame oe It 
My Spode a2 0427 1ty 


; -luwe : 
WSpod nN} = NV > "s fe m 
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Appendix D 


Model M-3 Optimized SWATH 
Table of Offsets 
and 


Design Worksheet 


by 


Model M-3 Table of Offsets 


Demihull Offsets 


TE ee Strut Offsets 


in inches 
shl2- ft. :2OFF(xh123 0° 32 
SiR 12 bodvl(x1)-ft 12 





R232) 


Se 


57) | Sa es 
-—~ } to 
2 - 
OV] 


role 
~~] 
pa) 
tJ 
<2) 


tJ 
w. 
tJ 
~] 
(St 


to 
L 
to 





1 re 
aK) 
~J | ft 


-_~ 
ox 
SS 


to] ro 
ry] 
—— 

=J 






tofro]to 
|| oe || ee | ee 
— |} — | WY Ua 
ere | ~~] a Le | 








to 





tJ 
=) 
ON 
tr 





to 
S 
re 
) 






ho 7 to 
ad ho 
OU tv 
Wf 
TH] Ww 









tol Nh 
sal Ihe 
CN 179G 
om mee) 
CN] to 


Le) 
WN 
ft. 
ho 


to] to to 
tO 
abla 
Oo] Oo] Cl— 
&e\| fe || = }| UA 





‘oD EP UN 
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Model M-3 The High Speed Swath 


This variant uses the same gross dimmensions as the baseline SLICE but uses a 
shape characteristic of a high speed swath. Tandem strut design was selected to 
minimize strut displacement. 


1a. Units Defined 


=e lh= 1M 


nm 25076 tt 


sec=1T Base Units 


Nautical Mile 


i ee Knot 
hr 
lton =2240 lb Long Ton 





1b. Fixed Variables 
2 =O 0662 
2 
Sec 
DP ayy = ae th a 
It 
i Regen eS it 
View oes HK ot ee 
sec 
9367 1] occa 
> Fy = 3 FO er 7 
{t 
ti 
| eset 
en SOS) = 
sec 


g=32 174+Itesec ? 


eee V2 0kt) Meoo kt 


Density cf Salt Water at S6° F An assumed average 
cperating temperature 


Kinematic Viscosity of Salt Water at coo F 


Density of Fresh Water at 68° F An assumed average tow 
tank termcerature. Use this if prediction is for a model test. 


Kinematic Viscosity of Fresn Water at 68° F 


Acceleration due to gravity 


Modify this to reflect the velocity range of interest. 
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1c. Indicator Settings 


Water 


$$ mW es “** indicates the snip is operating in Salt VVater 
"2" indicates the snip 1s operating in Fresh V/ater 
ae SAS TOC Ge 
09 =f SWel.p sw? fy; 
“a SW. : 
V alll Oe Dy Vi 
Body Segmentation 


Soe Seer 


2a. Hull Component and Required Displacement Definition 


as 177-Mon- 4” This defines the the reauireded (4) displacement of the vessel 


| , ae 
wy FAG igo The requireded (4) displaced volume 
pg 


2a.!. Strut Creation 


Strut Height 
Desired to allow demi hull submergence of at least one diameter. 


LY =X fla 
l] =6.355-e1n 
Strut Beam 
ve 6 By -325% Maximum Strut Beam Kequired Defind per the Baseline 
Considered a required parameter. 
B lex PSs on 


Overall Strut Length 
L, =241t. Overall Strut Length 


ere 


L , = 19.066 «in 
Length to Beam Ratio 


LB] =— LBI =7.385 Length to Maximum Beam Ratio 
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Length Breakdown 


aye pmb | = ely. Lenath of Parallel Midbody 
7 “Ee ate ae omb| Ly) =9.695-Strut Nose Length (2.4 x Dis optimum) 


L ro L \- Lomb | = ea Ly = ().695° Strut Tail Length (3.6 x D is optimum) 


Strut Nose and Tail Shape Exponent 


ere =e) OVER ee -7 ex 
nyo? Dips cr 


Calculation of Nose and Midbody Shépe 


\ | Se S22 le 








nt ~ = pmbi defines the x values of the subdivisions 
EN ns | 
<i 3 | alo defines the y values of the nose based 
a > oe ona parabola 
2 
‘omol( Sie ; ; 
San : defines the y values of the midsection 


ACN) Sat NE ary) sion te ne) 


Calculation of Tail Shape 


\1 =O ft seg Lay - Lomi - by defines the x values of the subdivisions 
; ; SAS 

Se oe Sage pmb . defines the y values based 

| Z on a parabola 


ae srt | 
4 


badvi(sli =at unl) Ll y= ee Alia ee Function that defines the 
l 


body geometry 


Striaceion 





2.a.ti. Wetted Surface and Volume Calculation For the Struts 


Perimeter 
Nose 
oi a ce mn 
3 > + | } 
NB ei ee mt 
ory BY 
Noy ee 
Tail Te: Lie Leen 16 LG 
Se a eae 7 t| | t | 
ee Ue, ery Sone EE ent), Yee 
| SLY By 
= a 
qT, ee leit 
Midsection 
Sp ==} amb! 
le a tt 
Total Perimeter 
= ia ee 
P , Ppt i 
P =3.24° 
Wetted Surface of the Strut 
WS. =PH = WS S51 71661 for one strut 
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Areas 


nese aled 


tail area 


mid section area 


Sq =B Jb omb| 


Total Area 


rea ie oe 


Volume per Strut 


sank = Area Fl 
Total Strut Volume 


VOL . =4 vol 


2.a.ti! Demi-Hull Creation 


- nae 
Ne alk 

Saene 
ee Onl 
ca ap: 
mea 


Area? =4) 242516 


- 
oO) = 2SeiG 


VOL , =0.513-i” 


Demi Hull Volume Required 





ue = lett for two hulls 


per demihull 


—— 
tr 
(J 


Dem! Hull Shape 


D h = 7 ft-2. Max diameter of the Hull D =o 3.36} mn 


Section #1-2 interface 


ad 


Dy] =f B Fe Oy ieee Or ae 


D i = 73-D . TD) h2 =4059e1n Section #2-3 interface 
[) h3 =) 1 D h3 =3.56] ¢in Section #3-4 interface 
1D hae = D 1 18 ha = 3.56] *1n Section #4-5 interface 
D> oe aa) : 2D ho O59 en Section #5-6 interface 
D ie, Weer Dia 2 a Section #6-/ interface 
Ee iia Overall Length of the null. eee? oT 
Li = (2.08 Ly Dt = 5.72¢mn Section #1 
ig 2 Soe iy ib 2? = 13 584 ein Section #2 
a ie eo) OS ie 1) Lk 3 ra 05 1) Section #3 
Lys eee is, ee = |7 874i Section #4 
Eas an) ee , is =) (eon Section #5 
h6 Soule i, I h6 =H 435-1n Section #6 
Lee, =a ale ae 1, Li = 13.584 en Section #7 
L. h : 
i 2D = 2357 
1) r 
Seament 1&2 
Calculation of Nose Cone Shape 
Deyn = es Nose cone shape exponent 


defines the x values of 


eee ot ee by 5 the subdivisions 


le 


ore ie \h1 2! mal 
vhi(shl2) -:=———. | - ———»_-__— 


> T 
e | 


Dp 1) 4 Se - 
OES Cs “hl . n2 hha; 
: - ly 








AlCl?) Sn ahi Beene it NIN as We aka 


Adding in Segment 3 


Ly) — Eh —Ey3, 


Divo D 143 ~ ike 











ca, Ea | 


=~ 
— —_ os nN 


Bl¢vhl2) 


Adding in Segment 4 


1112 at . 1 = ave 
Nnt2 UT seg Py by bag - bag, 


Das 
wns (sll?) ke 


+ 


Nhl2-L.. -=219 = seg; 


sat Nine SO etree I Rae es se 


the subdivisions 


defines the y values of the 
nose based on an ellipse 


~ 


defines the y values of 
- se2:,; section 2 


offsets for the first two body 
segments. 


defines the x values of the 
SUubDIVISIONS 
defines the y values 


of section 3 


offsets for the first three 
body segments. 


defines the x values of the 
subdivVISIONS 


defines the y values of 
section 4 


Ctx) ith Shy ee Bue — pees Page Ne 


Adding in Segment 5 


Shee eee 


~ 


<D Dassen 
vh5¢xhl2) - hf Shs 7 hd 








- 2 Lys 


re) 
ts 


Pye Overs Ete Gi le 


offsets for the first three body 
segments. 


defines the x values 
of the subdivisions 


MIZ~ Ly yp bya bag~bpg - ee. 


defines the y values of 
section 5 


Oi emer emai by ba yg, Cl(xh12),sn5(xh12)) 


offsets for the first three body 
segments. 


Adding in Seqment 6 


Nhl? =O0-it-seg.- Py) ~Lj5o-Ly3-ciy—-Ll,5-Ly6, defines the x values of the 


subdivisions 


~ 


ys ee nG 2 











5 : ~ | 
shots) = eee ae a eye Ey. ese |) 
- - Ling ; ” re os 
defines the y values of 
section 6 
TRIS ye 2: 


Ny foe ee = i = Ge - Liy4 _ Pew Sea Oe N12) 


offsets for the first three body 
ScCiems, 


Adding in segment 7 


| ae é defines the x values 
Shl2 SON see ee ee ae lee — Lie 


of the subdivisions 


GN Idee 





s ee Mat? ep Lp byg-bng- bps + lne Saco 
: ih 


defines the y values of 
section 7 


OPEC Sipe ee! Nf: a Bey _ Laie: - 3 _ ies > Lee ~ Pe ee Sn) sie SZ) 


offsets for the first three 
body segments. 


OPF( si 2) 





Wetted Surface of TWO demihulls 


Dales 
WS, -2- 2-4: OFF(xh12) dhl 2! WS , = 13.2+ft° 
- 0-11 
Volume of TWO demiuhulls 
< Ly 
VOLF = 2m x OFF(xh12)* dvhl2: VOL , = I.3+ft 
«0-11 , 
Volume Balance 
Available 
VOr Wali = me Ge) 1 - V¢ re s 
VOL Avani = 1 Sett 
Maya = Y OL avai Pg 
Navan! = | Ls S-Ib 
Required 
Vg=l sar 
excess Volume 
VOL ai =o (negative implies insufficient volume) 
Deck Area Lost compared to SLICE 
C San This is the combined ammount the bow/stern are 


cantelivered over the ends of the hulls. 


J 


ee we a> tee aes a — 
105 tt- —- C. 55-1 =660-t" (negative implies a gain in deck area) 
/. 


7) 


Model! Building Balance 
Bouyant Force 
ee Sas Fy, = 112.496-lb 
Weight Account 
FRAME 
W aoe |S Ib 


Force Block + Heave Post 


W hh] = aS ib 


Other Model Supplies 


W oth = 2-Ib 


Desired Reserve Bouyancy 


Allowed Weight of a gemihull + 2 Strut combo 


Fa- We- Wor WY oth > Ford 


W sp > 


Appendix E 
PrO|cet Pectenian 


and 


Test Journal 
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SLICE Project Test Plan for the Week of 6-10 March 1995 
Anticipated Number of Models to Test: 4 


Desired Data to Obtain: Resistance (x) for all 4 models, Dynamic Sinkage (z) for 2 
models. 


Model Descnptions: Model M-! is the LMSC baseline design. 

Model M-2 anticipates improvements to the LMSC baseline by 
improving the shape of the bodies but retaining the gross 
geometry of the baseline ship. 

Model M-3 1s an optimized SWATH based on the results 
presented by McGregor 

Model M-4 1s the Mid Foil Concept. 


Desired Testing Schedule: 


6 March 95 Set-Up M-1: M-1 Test 1 20runs 

7 March 95 M-1 Test 2; Set-Up M-2: M-2 Test | 30 runs 

8 March 95 M-2 Test l(finish) ; Set-Up M-3: M-3 Test 1  30runs 
9 March 95 M-3 Test 1 (finish); M-3 Test 2 30 runs 
10 March 9S Set-Up M-4; M-4 Test | 20 runs 


Test Description: 


Test 1: Resistance Locked in Heave and Pitch. 


Objective: To measure forces in the X and Z directions over the 
operating range of speeds. 


Discussion: Previous testing has indicated that a significant Z force is 
generated and that the magnitude of the force is geometry dependent. This method of 
testing replicates the LMSC methodology and therefore allows for a basis of model 
comparison with LMSC results. The LMSC justification for this methodology 1s that 
control surfaces will be employed to fix the operational draft. 


Description: With the model ngidly fixed to the carnage, the following 
sequence of speeds will be used to collect data for both X and Z forces : Model Scale 
Equivalent (MSE) of 3,6,9,12,16,20,25,30,35,40 kts for basic data. Each speed will 
accomplished with one run. Then MSE of 20,25,30,35,40 kts to demonstrate 
repeatability. If time in the day permits then additional runs will be made at MSE of 
10,11,13,14 and 22 kts to fill in hump data and the high speed knuckle. This 
represents a total of 20 runs for the test completion. 
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Test 2: Free to Heave. 


Objective: To quantify the difference in resistance results for a model 
rigidly fixed and one that 1s allowed to oscillate about a desired waterline. 


Discussion: Although control surfaces may hold the ship to a fixed 
depth in a gross sense of the word. it is unlikely that the system's sensitivity will 
exactly replicate a ngid mounting. This test will use a counterbalance to replicate the 
control surface but allow the model to oscillate about this point. 


Description: This test will use the Z force results from Test 1 to 
determine the counterbalance weight required. The model will be attached to the 
carnage free to heave but locked in all other respects. The necessary counterbalance 
weight be added and the model towed at the given speed. (Note, due to the small TPI 
associated with SWATH vessels, a method for adding the weight as the carriage 
comes to speed may be necessary. This is a Discussion Point for the 3 Feb meeting.) 
The test mg will require a positive stop in the negative z direction with a method of 
determining that the stop has been reached. This stop must be set to no more than 3 
inches above the design waterline of the model. ° 

The following sequence of speeds will be used to collect data for X 
force only with note being taken of the counterbalance weight added. Model Scale 
Equivalent (MSE) of 9,12.16,20.25,30,35,40 kts for basic data. Each speed will 
accomplished with one run. Then MSE of 20,25,30,35,40 kts to demonstrate 
repeatability. If time in the dav permits then additional runs will be made at MSE of 
3, 6,10,11,13.14 and 22 kts to fill in hump data and the high speed knuckle. This 
represents a total of 20 runs for the test completion. 
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index of Test Runs Conducted 6-10 March 1995 


Run# ModevVTest Full Scale Test Drag Lift 
Equivalent \ Velocity (Lf) (Lbf) 
(Kts) (FPS) 

1° M-1 in Baseline Config. 10 4.08 | 6.46 12a 
2: M-1 in Baseline Config. Sh Hi2t! 0.151 0.191 
M-1 in Baseline Config. 6} 2.41: 0.689 0.562 
M-1 in Baseline Config. 9. 3.64 : 4.24 1.09 
3: M-1 in Baseline Confia. 2 4.85 ; 9.08 2.425 
M-1 in Baseline Config. foo 5.47. 6.16 3.26 

6! M-1 in Baseline Config. 16 6.47 5.61 5.16 
7 M-1 in Baseline Config. 20 8.08 | 7.193 ee, 
R/7 M-1 in Baseline Config. 20 8.08 | FEA Te 
8| M-1 in Baseline Config. ha Void 10.66 10:53 
9: M-1 in Baseline Config. 30 V2 Von 19.54 15.52 
RQ M-1 in Baseline Config. 30 Zoe 19.39 15.46 
10' M-1 in Baseline Config. So 14.14 26.02 21.02 
11 M-1 in Baseline Config. 40 16.18. Gls GV 26.42 
12: M-1 in Baseline Config. 20) Ueto 621 1.02 
M-1 in Baseline Config. 11 4.45 8.91 1.64 

14) M-1 in Baseline Config. 16 S20. 6.86 2.74 
M-1 in Baseline Config. 14 5.66. 5/65 5.50 

16; M-1 in Baseline Config. 15 6.02 | 5.48 425 
17'_M-1 in Baseline Config. 40. 16.18 39.56 26.65 
18; M-1 with Spray Rail Aft 30 PAA 19.45 15.48 
19! M-4 with Spray Rail Aft 35 14.14 Va ino }e) 20.61 
20° _M-1 with Spray Rail Aft 40 16.18. 36.87 26.87 
24' M-2 in Baseline Config. 10 4.11 6.83 0.443 
25' M-2 in Basetine Config. an 1.22) OZ 0.247 
M-2 in Baseline Config. 6. 2.46 0.667 0.845 

27: M-2 in Baseline Config. 9. 3.68. 4.52 0.416 
M-2 in Baseline Config. 12 4.93 9.92 Zon 

28 M-2 in Baseline Config. 6 246 0.679 0.724 
29' M-2 in Baseline Config. mz 4.93 9.76 2.51 
30: M-2 in Basetine Config. 16 Gioc) 6.28 5.95 
31 M-2 in Baseline Config. 16) 6:58: 63 5.89 
32| M-2 in Baseline Config. 200 S235 TZ: 9.35 
R32 M-2 in Baseline Config. 2OP 8.23 | Woes 9.29 
34| M-2 in Baseline Config. 25 10.28. (8) 28 12.97 
R34 M-2 in Baseline Config. 20. 1078 10.19 13 00 
36| M-2 in Baseline Config. 30 AZiooe 14.54 V/.52 
R36 M-2 in Baseline Config. 30 12°35. VA°50 17.48 
38! M-2 in Baseline Config. Sle 14.4: 207 22/30 
R38 M-2 in Baseline Config. 35 14-4. 20.19 225 
40: M-2 in Baseline Config. 40 16.45 | 26.42 28.25 
R40 M-2 in Baseline Config. 40 16.45. 26.36 28.35 
42: M-2 in Baseline Config. 7 447 9.81 OL7 
43! M-2 in Baseline Config. qalde 472° 17S 1.82 
44° M-2 in Basetine Config. 1305 547. 6.68 B20 
45: M-2 in Baseline Config. eae 522 Thisl® 2.86 
46! M-2 in Baseline Config. 14.5 9.97 6.05 4.27 
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48 M-3 6! 258 0.801 1.01 | 
M-3 9! 3.88 4.22 1.98 | 
49 M-3 12 | 5.18. 6.08 1.08 | 
53 M-3 16! 6.91 8.24 4.48 | 
55 M-3 16 6.92 8.22 463. 
56 M-3 18) (als 8.93 6.4 | 
Sy, M-3 20) 8.65. 9.74 7.9 | 
58 M-3 Towed Backwards 75 3.48 3.17 1.26; 
59 M-3 25 10.81 13.87 Ze 
R59 M-3 Zoi, 10.81 13266 11.92: 
61 M-3 30} Ac 19.13 16.68. 
R61 M-3 30 | es 19.03 16.56 | 
63 M-3 Shey Weis 26 22.39 | 
R63 M-3 35 15:13 Zo15 22.19. 
66 M-3 40 17.29 So27 28.84 | 
R66 M-3 40 | 17.29 3001 28.7 
68° M-3 A485 4.97 eso eras 
69 M-3 14 5.97 6.74 227, 
70 M-1 in Baseline Config. 125 At 6.01 1.87 | 
M-1 in Baseline Config. 1125 458 7.14: 1.42. 
M-1 in Baseline Confa. 1k 478 G5: ee 
71 __M-1 in Baseline Config. 10° 4.05 6.37— 44138 
M-1 in Baseline Config. AS 468 9.43 2.29 
72 M-1 T+0.5" 1195 a7) 9.59 AS: 
M-1 T+1.0" 115 AT 1 
74 M-1 T+0.25" ies: 467 9.42 0.89 | 
R74 M-1 T+0.375" lien 47 9.45 0.4! 
76 M-1 Towed Backwards 20° 8.07 -9.91 10,22, 
7? M-1 Draft Adjustment 153 5.97 SeePlots See Plots . 
M-1 Draft Adjustment 15i 5.97 SeePiots See Plots 
78 M-1 Draft Adjustment ie) 5.97 SeePilots See Plots . 
M-1 Draft Adjustment 15) 5.97 SeePlots See Plots 
M-1 Draft Adjustment 15 5.97 SeePlots See Plots 
79‘ M-1 Draft Adjustment feo 6.97 SeePlots See Plots 
80 M-1 Draft Adjustment tivo 6.97 SeePlots See Plots 
81 —M-1 Draft Adjustment a TAle) 6.97 SeePiots See Pilots 
82 | M-1 In Line Config 4.05 O02. RATE 
M-1 In Line Config 4.85 9.07 2.67. 
83 M-1 In Line Config 1.21 0.158 Oo: 
M-1 In Line Config 2.41 0.693 0.6" 
84 M-1 In Line Config 3.64 4.41 1! 
M-1 In Line Config 6.47 5.74. 5.36. 
BS: M-1 In Line Config , 8.08 Lae 754 | 
86 M-1 in Line Config On 11.48 10.26 
87 | M-1 In Line Conf Peas, 22.29: i422) 
88 | M-1 In Line Config 14.01 SO. Ur 19.67 
89 M-1 In Line Config 16.19 44.51 26.45 
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Appendix G 


Data Scaling Worksheets 


This Calculation Scales the Test Data obtained for Model M-1 to Full Scale. 


NOTE: The System of Subscnpts in this Routine can become rather bulky... °R", "C", “M", "S" are Resistance, Coefficient 
of Resistance, Measured and Ship respectively... where measured refers to the measured model (Model M-2 Data 
unscaled) and Ship refers to the full scale calculation. Also, “s’, "p", "r’, “f, "t” as subscnpts are strut, pod, residual, 
frictional and total respectively. Other uses of symbois are defined as they anse. 


Input Data Section 


1:=0..16 Sets the data range 


RM a READ(mItlfx) Total measured X component force for test #1: fixed in pitch and heave. 
i 


VM. = READ(mIltlv) Speed increments used in the test. (inputted in fps) 
m= 17.26 Scaling Factor for Lengths 
. I me | | , 
VM kts. 7 VM.:-——— Measured Velocity in Knots (converts the inputs which are in fps) 
i ' 1.688 
= 
Weds = YMigg Va Scaled Velocity in Knots 
1 1 
i 
VS. = VM.-jA Equivalent Scaled Velocity in FPS 
a2 
Ure se 1.2260-10 R*2/sec Kinematic Viscosity of Fresh Water at 59 F (Test Temp) 
= -5 ; . 
Vow = 1.1794-10 fi*2/sec Kinematic Viscosity of Salt Water at 65 F (Assumed Operating Temp) 
Dia =. 2004 Ib sec*2/ft*4 Density of Fresh Water at 59 F (Test Temp) 
nS 1.9890 Ib sec*2/ft*4 Density of Salt Water at 65 F (Assumed Operating Temp) 


Geometry Section 


LM, =1.39 
SM, =0.715 
LM, =2.14 
SM, = 2.01 


SM, =4:SM,+4-SM,, 


Resistance Calculations 


feet Length of one Strut 


R42 Wetted surface of one Strut 


feet Length of one Pod 


ft*2 Wetted Surface of one Pod 


ft*2 Wetted surface of the entire underwater body 


NOTE: In this section Rn is calculated for the pod and the strut separately since their lengths 
are different. Coefficients of friction are then calculated separately for each based on the 1957 
ITTC Line. These are converted to Resistance Values and then subtracted from the Total 
Resistance value. Finally this is converted to a Coefficient of residuary resistance which is equal 


for the model and the ship. 


VM.-LM 
RnM 5 = eet ae es 
i v Ay 
VM..-LM 0 
RnM p. = 
i U fy 
0.075 
OO = : 
floe(RaM : ~ 2| 
1 / 
CM : 0.075 


operating Rn of the Strut 


operating Rn of the Pod 


Coefficient of Fnctional Resistance based on Strut Length 
and using the 1957 ITTC line. 


Coefficient of Fnctional Resistance based on Pod Length 
and using the 1957 ITTC line. 
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Frictional Resistance of four Struts 


(VM,)*-4:SM , 


RM ¢, = CM ¢5 -—-2 fy 


Frictional Resistance of four Pods 


-P fw (VM,)"4°SM | 


fs 


Residual Resistance for the Model as measured. 


& ——— 
3 = 
1 a 
2 gis 
z ig) 








Note: these values for Rn are Low for ITTC 1957 use. Should be on order of 10*7 
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Scaled Resistance Calculations 


CS 5 = CM : Residual Resistance Coefficient remains the same in the saling. 
i i 
VS-LM <3 
1 : 
RnS S = ————____—— operating Rn of the Strut in the new scale 
Y sw 
VS.-LM ..-4 
= Se) P . . 
Sa operating Rn of the in the new scale 
Rns |, ting Rn of the Pod in th | 
“Sw 
0.075 | | | 
GS {S. 332 Coefficient of Fictional Resistance based on scaled Strut Length 
' slog; RnS ,\ - 2\? and using the 1957 ITTC line. 
| 7 
| 
0,075 | _ . 
Gs fp = Coefficient of Fnctional Resistance based on scaled Pod Length 
. ( log/ RnS _ \- ne and using the 1957 ITTC line. 
Da 
CS an 0.0004 Correlation Allowance 
Roe Coe = (VS ee SM ae Frictional Resistance of the Struts 
ts, fs. Paw 2 


RS fp. =U fp? SW (VS.)*4-SM wie Fnctional Resistance of the Pods 


AVIA . ae . 
P sw (vs,\* SM pA Scaled Residual Resistance 


RS a cS a sw (VS.)"SM ee Correlation 


. 


RS ; SRS ts) SS pie spate 


Scaled Scaled 
Speed Speed 
(fps) (knots) 


WRITE(mlfulfx) =RS y 


RS fi ZS fp. iS ts. Pe 


i 


RS 


a. 


CA. = 
© RS le 
i 


WRITE(mfulfp) = FF. 


44860.57 


163618.18 
























488.89 
680.46 
0480.51 
0929.77 
1798.85 
4370.38 
6882.96 
41361.02 


18667.49 


31883.98 


1 1 







27 Nod 
8777.95 
5868.95 
4672.03 
0738.84 
$766.63 
8587.27 
33802.96 





PodF rictional 


Total M-2 Residual 
Resistance Component Component 
at Full 
Scale 
RS, RS, 


RS gy 





WRITE(mIfulv) = VS yy. 


WRITE(m I fulrp) = RF. 


StrutFrictional 
Component 


RS fs. 





Reynolds # 
of the Strut 


























Reynolds # 
of the Pod 

Ra , RnS 5, 

15743425 
31356738 
47360385 
52694933 
53085266 
5789937 | 
60891923 
63103809 
68308247 
71170688 
73642796 
5467882 | 84181783 
10512964 
13154216 
15782457 
18397687 
21051951 





Writes the scaled data to a file for future use. 


percent of total resistance that is attributable to correlation allowance 


WRITE(mIfulcap) = CA, 


Percent of total resistance 
attributable to friction/wave. 


This Calculation Scales the Test Data obtained for Model M-2 to the scale that Model M-1 
is built to. This enables a direct comparison of the Data obtained from the two models. 


NOTE: The System of Subscnpts in this Routine can become rather bulky... "R*, "C*, °M’, "S" are Resistance, Coefficient 
of Resistance, Measured and Scaled respectively... where measured refers to the measured model (Model M-2 Data 
unscaled) and Scaled refers to the Model M-2 data scaled to the Model M-1 scale factor. Also, “s", "p", °r, “f, "t* as 
subscripts are strut, pod, residual, frictional and total respectively. Other uses of symbols are defined as they arise. 


Input Data Section 


RM, =READ(m2tIfx) 


VM. =READ(m2tlv) 


, _ 16.68 
17.26 
VM iis = VM:—— 


' 1.688 


VS kts. = Mi kts. - /, 
ie 
VS, = VM.vh 


D Ry = 1.2260-10° 


= -5 
Deel 40 
P Ay = 1.9384 
Ope a 1 ed 


Sets the data range 


Total measured X component force for test #1: fixed in pitch and heave. 


Speed increments used in the test. (inputted in fps) 


Scaling Factor for Lengths M2/M1 scale factor. 


Measured Velocity in Knots (converts the inputs which are in fps) 


Scaled Velocity in Knots 


Equivalent Scaled Velocity in FPS 


ft*2/sec Kinematic Viscosity of Fresh Water at 59 F (Test Temp) 


ft*2/sec Kinematic Viscosity of Salt Water at 65 F (Assumed Operating Temp) 


lb sec*2/ft*4 Density of Fresh Water at 59 F (Test Temp) 


lb sec*2/ft*4 Density of Salt Water at 65 F (Assumed Operating Temp) 


Geometry Section 


ON RYE 
SM , = 1.059 
LM , =2.32 

SM , = 2.642 


SM, =4:SM,-4-SM, 


Resistance Calculations 


feet Length of one Strut 


ft*2 Wetted surface of one Strut 


feet Length of one Pod 


ft*2 Wetted Surface of one Pod 


ft*2 Wetted surface of the entire underwater body 


NOTE: In this section Rn is calculated for the pod and the strut separately since their lengths 
are different. Coefficients of friction are then calculated separately for each based on the 1957 
ITTC Line. These are converted to Resistance Values and then subtracted from the Total 
Resistance value. Finally this is converted to a Coefficient of residuary resistance which is equal 


for the model and the ship. 


VM--LM . 
iS 
‘i ive 
VM;-LM 5 
RnM, =—— 
P; VD Ay 
0.075 
CM 5, = as «ee 
(eee | 2 
CM oe 0.075 
(log(RnM p)- 2° 


operating Rn of the Strut 


operating Rn of the Pod 


Coefficient of Frictional Resistance based on Strut Length 
and using the 1957 ITTC line. 


Coefficient of Frictional Resistance based on Pod Length 
and using the 1957 ITTC line. 


Frictional Resistance of four Struts 


VM,)*-4-SM . 
V/ 


( 


| 


‘= 70 . 
3 fw 


CM fe 


Frictional Resistance of four Pods 


45M 5 


(VM, )* 


PAW 


Residual Resistance for the Model as measured. 


- RM, 
SM, 


1 
26 


P fw | 


RM 1). - RM gg, 
RM , 
VM. } 


RM , 
eis = 














5 


Note: these values for Rn are Low for ITTC 1957 use. Should be on order of 10*7 


Scaled Resistance Calculations 


Residual Resistance Coefficient remains the same in the saling. 
1 i 


VS-LM oh 
= J 


Rns . operating Rn of the Strut in the new scale 
VS.-LM p A 
RnS p ————— operating Rn of the Pod in the new scale 
; Paine 
. 0.075 i” 
Cx fs = Coefficient of Frictional Resistance based on scaled Strut Length 
i ‘log/RnS 2.5 ae and using the 1957 ITTC line. 
1: 
0.075 tj | 
Cy fa Coefficient of Frictional Resistance based on scaled Pod Length 
‘log/RnS p.7 2: 2 and using the 1957 ITTC line. 
\ ; i 
Cs) Correlation Allowance 
a 
RS =CS ¢ an ew ie Frictional Resistance of the Struts 
fs. ; ts. 5 fw ae sae 


RS fp. =CS be Ay VS,)°:4:SM me Frictional Resistance of the Pods 


oe. 2 2 
Soe ry Pw VS\"SM yA 


Scaled Residual Resistance 


RS t. 


“RS g. + RS gy -RS_ 


I 
— 
tI 


Scaled Scaled Total M-2 Residual PodFrictional StrutFrictional Reynolds # Reynolds # 


Speed Speed Resistance Component Component Component of the Strut of the Pod 
(fps) (knots) at 17.26:1 
VS VS ee RS : RS fp. 


0.09 
0.91 
ae 
1.08! 
1.81) 
7.38 
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YJ | — 
LJ | — 
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WRITE(m2sclfx) =RS t. WRITE(m2sclv) = VS. Writes the scaled data to a file for future use. 
Sos 
RS aR peek FE. == — RE 22.) 2.5 Percent of total resistance 
ft. fp. fs. 1 ' 1 1 ae 
‘ : ! RS t attributable to friction/wave. 


WRITE(m2ff) He WRITE( m2rt) =RF. 


tO 
—=— 
ws) 


This Calculation Scales the Test Data obtained for Model M-2 to Full Scale. 


NOTE: The System of Subscripts in this Routine can become rather bulky... "R"”, °C", "M", "S° are Resistance, Coefficient 
of Resistance, Measured and Ship respectively... where measured refers to the measured modei (Mode! M-2 Data 


unscaied) and Ship refers to the full scale caiculation. Also, “s", "p", °r, "f, "t’ as subscripts are strut, pod, residual, 
frictional and total respectively. Other uses of symbols are defined as they arise. 


Input Data Section 


i-0,, 15 Sets the data range 


RM t] = READ(m2tlfx) Total measured X component force for test #1: fixed in pitch and heave. 
i 


VM. = READ(m2t i) Speed increments used in the test. (inputted in fps) 
4 = 16.68 Sealing Factor for Lengtns M2/M1 scale factor. 
_ | , 
VM kts. ~ ON ras Measured Veiocity in Knots (converts the inputs which are in fps) 
i 1.688 
VS kts. ~ VM kts.V A Scaied Velocity in Knots 
1 i 
VS, = VM. A Equivalent Scaled Velocity in FPS 
5 
Vay = }.2260-10 fR*2/sec Kinematic Viscosity of Fresh Water at 59 F (Test Temp) 
-5 
Vow = 1.1794-10 ARA2/sec Kinematic Viscosity of Salt Water at 65 F (Assumed Operating Temp) 
P fy = 1.9384 Ib sec*2/ft*4 Density of Fresh Water at 59 F (Test Temp) 
Pie - 1.2890 Ib sec*2/ft*4 Density of Salt Water at 65 F (Assumed Operating Temp) 


Geometry Section 


LM , = 1.4375 feet Length of one Strut 

SIV era aOoe ft*2 Wetted surface of one Strut 
LM p= 2.32 feet Length of one Pod 

SM p = 2.642 ft*2 Wetted Surface of one Pod 


SM fos 4-SM Agi 4-SM p ft*2 Wetted surface of the entire underwater body 


Resistance Calculations 


NOTE: In this section Rn is calculated for the pod and the strut separately since their lengths 
are different. Coefficients of friction are then calculated separately for each based on the 1957 
ITTC Line. These are converted to Resistance Values and then subtracted from the Total 


Resistance value. Finally this is converted to a Coefficient of residuary resistance which is equal 
for the model and the ship. 


VM.-LM , 
RnM S. 7 operating Rn of the Strut 
tay 
VM. LM p 
RnM p. = ——— operating Rn of the Pod 
a i 
0.075 
CM fo. = Coefficient of Frictional Resistance based on Strut Length 
: (log RnM 5 Ve 9\2 and using the 1957 ITTC line. 
\ SH 
0.075 
CM oe Coefficient of Fictional Resistance based on Pod Length 
: fara p ay, \2 and using the 1957 ITTC line. 
\ 
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Frictional Resistance of four Struts 


4-SM , 


=e iy VM, )* 


1 


| 


CM ¢ 


RM fs. 


Frictional Resistance of four Pods 


(VM,)*4-SM ,, 


—"p ry 
9 fw 


I 


= CM f, 


BOs 
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Residual Resistance for the Model as measured. 


—RMp 
2 
SM ; 


1 
P Aw’ | 


i 


RM, =RMy) -RMg 
i 
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Note: these values for Rn are Low for 1TTC 1957 use. Should be on order of 10°7 


Scaled Resistance Calculations 


CS - = CM “ Residual Resistance Coefficient remains the same in the saling. 
1 1 
VS--LM oA 
1 é 
RnS eer aaa operating Rn of the Strut in the new scale 
t 18) SW 
VS.LM 74 
i Pp 
RnS p -— operating Rn of the Pod in the new scale 
l Uv Vy 
0.075 | 7 | 
CS fs a Coefficient ofFrictional Resistance based on scaled Strut Length 
! log) RnS _ ae and using the 1957 !TTC line. 
i] / 
0.075 
CS {G.. ° oS a ee Coefficient of Frictional Resistance based on scaled Pod Length 
log; RnS p _ ae and using the 1957 ITTC line. 
: , 
CS , = 0.0004 Correlation Allowance 
RS =CS. cn ey SM He Frictional Resistance of the Struts 
fs. ~~ fs.5 Pow ey 8 aes 


— 


RS fp. = 08 pas Pace (VS,)\*-4. Si) “rictionalResistanesen neleees 


/ P 
Csr VS.)?.SM pA? 7 7 
RS r =CS 15 Psw'( : | SM yA Scaled Residual Resistance 
=CS.- VS.)?.SM ph 7 
Ko a. a len a P sw’ ; m. tA Correlation 


a. 
1 


RS, =RS pg, —RSg + RS, - RS 


Scaled Scaled Total M-2 


Speed Speed Resistance 
(fps) (knots) at Full 
Scale 


WRITE(m2fulfx) =RS, 


RS ¢ eRS ipo RS ts. 


i 


RS 


a. 
fo 
RS _. 





WRITE(m2fulfp) =FF. WRITE(m2fulrp) = RF. 





WRITE(m2fulv) = VS pH. 


PodFnctional StrutFnctional 


87.1] 
14.68 
Bike 
09.02 
44.46 
044.18 
bist 
25739 
370.81 


7 


200.31 | 
727.02 


025.51 
29173 
4404.6 

189.62 

231.34 
10540.35 
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percent of total resistance that is attributable to correlation allowance 










10129794 
042565 | 
0555445 
4125783 
7114902 
9190680 


43342235 
5418013 


4634465 
8334597 
SS a 1: 
102543410 
1956478 
136586162 
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WRITE( m2fulcap ) = CA. 


Reynolds # Reynolds # 
Component Component of the Strut of the Pod 
RS fp, RS fs. RnS s RnS p 
















16348607 
2965225 
Pos 8oe 
5076046 
9900225 
3250349 
6064454 
9950598 
3300723 
0000972 
8175275 
110286097 
SSI 
165496147 
192207 1G 
220438188 
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Writes the scaled data to a file for future use. 


Percent of total resistance 
attributable to friction/wave. 


EHP ship 


oa 


VS. 


M-2 in full scale 





1049.079 








This Calculation Scales the Test Data obtained for Model M-3 to the scale that Model M-1 
is built to. This enables a direct comparison of the Data obtained from the two models. 


NOTE: The System of Subscripts in this Routine can become rather bulky... "R®, °C°, "M", "S" are Resistance, Coefficient 
of Resistance, Measured and Scaled respectively... where measured refers to the measured model (Model M-3 Data 
unscaled) and Scaled refers to the Model M-3 data scaled to the Model M-1 scale factor. Also, “s*,"p", "rr, "f, "t" as 
subscnpts are strut, pod, residual, frictional and total respectively. Other uses of symbols are defined as they arise. 


Input Data Section 


p= 0..16 Sets the data range 


RM os READ( m3t ] fx) Total measured X component force for test #1: fixed in pitch and heave. 
1 


VM. =READ(m3tiv) Speed increments used in the test. (inputted in fps) 





Scaling Factor for Lengths M3/M1 scale factor. 


i726 
| ace | 
VM kts. 7 VM.:-——— Measured Velocity in Knots (converts the inputs which are in fps) 
' 1.688 
a 
VS kts. ~ VM kts. V ). Scaled Velocity in Knots 
1 | 
VS, = VM.-4) A Equivalent Scaled Velocity in FPS 
‘5 
Vey = 1.2260: 10 ft*2/sec Kinematic Viscosmy of Fresh Water at 59 F (Test Temp) 
-§ 
ays 1.1794-10 ft*2/sec Kinematic Viscosity of Salt Water at 65 F (Assumed Operating Temp) 
Pay = 1.9384 Ib sec*2/ft*4 Density of Fresh Water at 59 F (Test Temp) 
Pow = 1.9890 Ib sec*2/ft*4 Density of Salt Water at 65 F (Assumed Operating Temp) 


t 
tJ 
= 


Geometry Section 


LM s 7 1.589 feet Length of one Strut 

SM 5 Las ft*2 Wetted surface of one Strut 
LM p = 35,955 feet Length of one Demihull 

SM p> 5.45 ft*2 Wetted Surface of one Demihull 


SM Ci 4+.SM st 2-SM p ft*2 Wetted surface of the entire underwater body 


Resistance Calculations 


NOTE: In this section Rn is calculated for the pod and the strut separately since their lengths 
are different. Coefficients of friction are then calculated separately for each based on the 1957 
ITTC Line. These are converted to Resistance Values and then subtracted from the Total 


Resistance value. Finally this is converted to a Coefficient of residuary resistance which is equal 
for the model and the ship. 


VM--LM , 
RnM i operating Rn of the Strut 
Painy 
VM. LM Dp 
RnM p =... operating Rn of the Demihull 
1 vO fw 
‘ 0.075 
CM fy Se Coefficient of Frictional Resistance based on Strut Length 
log/RnM z Ve he and using the 1957 ITTC line. 
i} 
0.075 
CM fp ae Coefficient of Frictional Resistance based on Demihull Length 
i log/RnM ' : oh and using the 1957 ITTC line. 
. i / 


ty 
ty 
—— 


Frictional Resistance of four Struts 


VM, )*:4-SM , 


l 
RM ¢, = CMs. =p py 


Frictional Resistance of two Demihulls 


VM, )*:2:SM ,, 


fw 


| 


aap 


Residual Resistance for the Model as measured. 


-RM¢p 
SM jf} 


i 
4 


fy" 


RM; - RMg 
RM, 
VM 


RM, 
Chee 


i 





Note: these values for Rn are Low for ITTC 1957 use. Should be on order of 10*7 


™N 
Cot 


Scaled Resistance Calculations 


Cs 2 (Ohl Residual Resistance Coefficient remains the same in the saling. 
Yr. r 
1 i 
VS-LM .2 
RnS = operating Rn of the Strut in the new scale 
S 
VS. LM p d 
RnS ° = ——_____-— operating Rn of the Pod in the new scale 
: 0.075 _ | 
CS fs a Coefficient of Frictional Resistance based on scaled Strut Length 
' log) Rn . 5 4 and using the 1957 ITTC line. 
| . 
O07> . a . 
> coefficient of Frictional Resistance based on scaled Pod Leng 
CS = Coefficient of Frictional Resist based led Pod Length 
log) RnS p ae : and using the 1957 ITTC line. 
Coa Correlation Allowance 


ee. cone ee 
RS fs. =CS isa fx" | MS SM ok Frictional Resistance of the Struts 


-_ 


RS fp. =CS 5 fw (VS,)*.2-SM wo Frictional Resistance of the Demi Hulls 


Roar ses 


. \ q _ 2 
Pau. (VS. MA a Scaled Residual Resistance 
i i 


to | — 


RS, =RS p+ RS gy +RS, 


1 


tI 
to 
Lad 


Reynolds # 
of the Pod 


Reynolds # 
of the Strut 


PodFnetional StrutFrictional 
Component 


Component 


Residual 


Total M-3 


Component 


Resistance 


at 17.26:1 
Scale 
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Percent of total resistance 
attributable to friction/wave. 


Writes the scaled data to a file for future use. 
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WRITE( m3rf) 


EF. 


WRITE( m3 ff) 


ie | 


This Calculation Scales the Test Data obtained for Model M-3 to full scale. 


NOTE: The System of Subscripts in this Routine can become rather bulky... "R","C", "M", "S" are Resistance, Coefficient 
of Resistance, Measured and Scaled respectively... where measured refers to the measured model (Model M-3 Data 
unscaled) and Scaled refers to the Model M-3 data scaled to the Model M-1 scale factor. Also, "s","p", "r, "f, "t" as 
subscripts are strut, pod, residual, frictional and total respectively. Other uses of symbols are defined as they arise. 


Input Data Section 


RM, = READ(m3tlfx) 


VM. =KREADC mtv) 


Y Ay 


SW 


P fw 


P sw 


1 


2c 


1.1794-10°° 


1.9384 


1.9890 


Sets the data range 


Total measured X component force for test #1: fixed in pitch and heave. 


Speed increments used inthe test. (inputted in fps) 


Scaling Factor for Lengths M3/M1 scale factor. 


Measured Velocity in Knots (converts the inputs which are in fps) 


Scaled Velocity in Knots 


Equivalent Scaled Velocity in FPS 


ft*2/sec Kinematic Viscosity of Fresh Water at 59 F (Test Temp) 


ft*2/sec Kinematic Viscosity of Salt Water at 65 F (Assumed Operating Temp) 


Ib sec*2/ft*4 Density of Fresh Water at 59 F (Test Temp) 


lb sec“2/ft*4 Density of Salt Water at 65 F (Assumed Operating Temp) 


td 
ho 
Ln 


Geometry Section 


LM , = 1.589 
SM , = 1.716 
LM, =e 
SM 16115) 
SM, =4SM,-2:SM, 


Resistance Calculations 


feet Length of one Strut 


ft*2 Wetted surface of one Strut 


feet Length of one Demihull 


ft*2 Wetted Surface of one Demihull 


R*2 Wetted surface of the entire underwater body 


NOTE: In this section Rn is calculated for the pod and the strut separately since their lengths 
are different. Coefficients of friction are then calculated separately for each based on the 1957 
ITTC Line. These are converted to Resistance Values and then subtracted from the Total 
Resistance value. Finally this is converted to a Coefficient of residuary resistance which is equal 


for the model and the ship. 


VM.-LM , 
RnM. =————— 
| ais 
VM.:-LM p 
Rv =e 
- » fw 
0.075 
CM fs = 5 
‘log/RnM . | - 2 
O:075 
CM fp. = ; 
‘log/RnM p,) 7 2 


operating Rn of the Strut 


operating Rn of the Demihull 


Coefficient of Frictional Resistance based on Strut Length 
and using the 1957 ITTC line. 


Coefficient of Frictional Resistance based on Demihull Length 
and using the 1957 ITTC line. 


Frictional Resistance of four Struts 


VM. \*:4-SM , 


| 


peal ny) 


RM ¢, = CM ¢, 


1 


Frictional Resistance of two Demihulls 


. 25, : 
VM; °2-SM , 


ee ian 


1 


| 


Residual Resistance for the Model as measured. 


—RMs, 
2 
SM , 








Note: these values for Rn are Low for ITTC 1957 use. Should be on order of 10*7 


Scaled Resistance Calculations 


CS as CM . Residual Resistance Coefficient remains the same in the saling. 
1 1 
VS-LM oh 
| 1 
RnS a operating Rn of the Strut in the new scale 
; U sw 
VSO 2 
ea Pp | 
RnS oo operating Rn of the Pod in the new scale 
U sw 
ae OO7> . . | 
cs — = Coefficient of Frictional Resistance based on scaled Strut Length 
' ‘log! RnS 517 me and using the 1957 ITTC line. 
\ | 


1/ 


: 0.075 | | 
coefficient of Fictional Resistance based on scaled Pod Leng 
CS fy Coefficient of Frictional Resist based led Pod Length 
: ‘log; RnS ang Ne and using the 1957 ITTC line. 
i/ / 
GS Ae 0.0004 Correlation Allowance 


a 2a cng 32 
RS fs. = cS ina a VS.) -4-SM g * Frictional Resistance of the Struts 
1 1 


RS fp, =CS fp? ae AS. \° 2-5M ae Frictional Resistance of the Demi Hulls 


P 
| roe wae ; 
Rs , =CS ay P oy’ (VS.)*.SM t = Scaled Residual Resistance 
Romesh Vs 5 OM Correlati 
a. : a, Psw ( S;) OVE A orrelation 


RS, =RS pg. ~RS gy ~ RS, +RSy 


tO 
I 
00 


Scaled Scaled Total M-3 

Speed Speed Resistance 

(fps) (knots) at Full 
Scale 


VS kts. 





WRITE( m3fulfx) =RS " 


os ft. =RS fp. ae RS fs. 





Residual 
Component 


WRITE(m3fulv) = VS p< 


PodFrictional 
Component 


RS fp 

























StrutFrictional Reynolds # Reynolds # 

Component of the Strut of the Pod 
RS fs Ris . Ras p 
76475883 
1031535] 
11501024 
12360631 
36207094 13575951 
[4168787 
14731982 
15354460 
17696163 
20497315 
73061332 
25640170 
39493606 
5854.09 32042802 
38534359 
119609899| [44848066 


Writes the scaled data to a file for future use. 





Percent of total resistance 
attnbutable to friction/wave. 


percent of total resistance that is attributable to correlation allowance 


WRITE(m3fulcap) = CA. 


= RF. 


WRITE( m3fulrp ) 


WRITE(m3fulfp) = FF. 


M-3 in full scale 
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Appendix H 


Approximation of Impingement Force 


This Worksheet approximates the impingement force created by 
the rooster tail on the towing frame. 


64 
p 3 


V ve = 10-— 
sec 


ee Sec 


ane 40° ft sec 


vol = 1.25-in-6-in- (Vt 
vol = 0.13+ft° 
Mm, =p-vol 
Meo oop 
F impact ~™Mw'4x 
F impact = 10:36-lbf 


density of water 


estimate of the maximum relative x component 
of velocity between the frame and water before impact. 


water x component of velocity after the impact. 


an approximaton for the time of impact 
of the water on the frame 


acceleration of the water 


an approximation for the volume of water that Is 
impacted during the impact period. 


mass of the water that is impacted dunng the impact 
penod. 


Approximated Impact Force on the frame. 
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Appendix. | 


Numerical Codes 


28s 


FKKKKKKEKES Oana 


% This routine solves for the total resistance of an 8 bodied SLICE tvpe 
% vessel. Solution to the friction resistance is based on the ITTC 1957 

% friction line. Solution to the wave resistance portion is based on the 

% solution to Mitchell's Integral proposed by Lunde in 1951 and expanded 
% by Wilson. 


% GENERAL INPUTS SECTION 
Is= input(‘Input the slowest integer speed as x.xx (f/s) ="): 


hs=input(‘The number of increments='): 
lani=input(‘The incrementation (decimal)='); 


di | 


water=input('For Fresh Water enter |. Salt Water enter 2: '); 
if water==1: 


rho= 1.9367: 
nu= 1.0804E-5: 
elseif water==2: 
rho= 1.9912: 
nu= 1.3343E-5: 
else 
disp(’ ') 
disp(’ ') 


Gite tS tte eet eS e ee Tee T LEEKS SS See eee es!) 


disp(‘A valid entry was not made for water type’) 


disp(The program may not provide valid results’) 
i i hn eager x oe ee EEE TESS Ee) 


end 
end 
disp(’ ') 
disp(’ ’) 
disp(’ ') 
disp('‘The bodv numbering 1s: |=fwd stbd strut’) 
disp(' 2=fwd port strut’) 


234 


disp(' 3=aft stbd strut’) 

disp(’ 4=aft port strut’) 

disp(’ 5=fwd stbd pod’) 

disp(’ 6=fwd port pod’) 

disp(’ 7=aft stbd pod’) 

disp(’ 8=aft port pod’) 

disp( ‘) 

disp( ’) 

disp( ’) 

disp(‘The global coordinate svstem 1s:') 

disp( X positive in the direction of motion’) 
disp(' VY positive to port’) 

disp(‘ Z positive upward ) 

disp(’ The origin 1s at the undisturbed free surface’) 


TS= input(‘The height of the struts (ft)='): 
LS= input(‘The length of the struts (ft)='): 
°oWSS=input(‘The wetted Surface of one strut (ft®2)='); 


cs6=input(‘The 6th power coefficient for the strut potential fn='): 
ces5=input(The 5th power coefficient for the strut potential fn='): 
cs4=input(‘The 4th power coefficient for the strut potential fn='): 
cs3=input(‘The 3rd power coefficient for the strut potential fn='): 
cs2=input( The 2nd power coefficient for the strut potential fn='): 
csl=input('The Ist power coefficient for the strut potential fn='); 
csO=input(The Oth power coefficient for the strut potential fn=); 


TP= input(‘The max diameter of the pods (ft)="); 
LP= input(‘The length of the pods (ft)='): 
°%WSP=input('The wetted surface of one pod (ft*2)="): 


cp6=input('The 6th power coefficient for the pod potential fn=’): 
cp5=input(‘The 5th power coefficient for the pod potential fn=’): 
cp4=input('The 4th power coefficient for the pod potential fn='); 
cp3=input(‘The 3rd power coefficient for the pod potential fn='); 
cp2=input('The 2nd power coefficient for the pod potential fn='); 
cpl=input(‘The Ist power coefficient for the pod potential fn='); 
cp0=input(‘The Oth power coefficient for the pod potential fn=’): 


sb]=input(‘The x setback of body | (+ or - ft)="'); 
osl=input(‘The y offset of body | (+ or - ft)=): 
subl=input(‘The z submergence of body! (— or - ft)="): 


tO 
(2 
ta 


sb2=input(‘The x setback of bodv 2 (- or - ft)='): 
os2=input('The v offset of body 2 (— or - ft)="'): 
sub2=input('The z submergence of bodv2 (+ or - ft)='): 


sb3=input('The x setback of body 3 (- or - ft)='): 
os3=input(‘The v offset of bodv 3 (- or - ft)='); 
sub3=input("The z submergence of bodv3 (+ or - ft)='): 


sb4=input(‘The x setback of body 4 (- or - ft)='): 
os4=input(‘The y offset of body 4 (— or - ft)='): 
sub4=input(‘The z submergence of body4 (+ or - ft)="): 


sb5=input(‘The x setback of body 5 (- or - ft)="): 
oss=input('The v offset of body 5 (— or - ft)='): 
sub5=input(‘The z submergence of bodv5 (+ or - ft)=): 


sb6=input(‘The x setback of body 6 (— or - ft)='): 
os6=input('The v offset of body 6 (— or - ft)='): 
sub6=input(‘The z submergence of bodv6 (+ or - ft)="): 


sb/=input(‘The x setback of body 7 (- or - ft)='): 
os/=input('The v offset of body 7 (— or - ft)='): 
sub7=input('The z submergence of body? (+ or - ft)="): 


sb8=input('The x setback of body 8 (- or - ft)='): 
os8=input(‘The v offset of body 8 (+ or - ft)="). 
sub8=input(‘The z submergence of bodv8 (+ or - ft)='): 


tseg=input("The number of Theta Segments='); 
zseg= input(‘The number of segments in Z='); 
xseg= input(‘The number of segments in X=’); 
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*o This section calls the routines that computes the Rw value 


° for each of the body and interface terms and the Rf value for 
* the body terms only. 


Rwb=0: 

Rw11=0: 
Rw22=0: 
Rw33=0: 
Rw44=0. 
Rw55=0: 
Rw66=0: 
Rw/7=0: 
Rw88=0: 
Rws=0: 

Rw15=0: 
Rw26=0: 
Reon 
Rw48=0: 
Rwi=0: 

Rw13=0: 
Rv )7=0. 
Rw53=0: 
Rw57=0: 
Rw24=0: 
Rw28=0: 
Rw64=0: 
Rw68=0: 
Rwd=0: 

Rw14=0: 
Rw 1 8=0: 
Rw54=0: 
Rw58=0: 
Rw23=0: 
Rw27=0: 
Rw63=0: 
Rw67=0: 
Rwt=0: 

Rw12=0: 
Rw 16=0: 
Rw52=0: 
Rw56=0: 
Rw34=0: 
Rw38=0: 
Rw74=0: 
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Rw78=0: 
Rf=0; 
Rfs=0; 
Rfp=0: 
Rt=0: 


OE a aia BODY RERMWISt ttt ++ tt + 
oneone 

twotwo 

thrthr 

forfor 

fivfiv 

SIXS1X 

sevsev 

ateate 


ie SAME CORNER INTERACTION TERMS****** 
onefiv 
tWOS1X 
thrsev 
forate 


Pee LONGITUDINAL INTERACTION TERMS********* 
onethr 

onesev 

fivthr 

fivsev 

twofor 

twoate 

six for 

sixate 


Paes DIAGONAL INTERACTION TERMS********* 
onefor 

Ooneate 

fivfor 

fivate 

twothr 

twosev 

sixthr 

SIXSEV 


pee elie No Ven SE DIAGONAL TERMS*********** 
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onetwo 
onesix 
fivtwo 
fivsix 
thrfor 
thrate 
sevtfor 
sevate 


Rwb=Rw11+Rw22-Rw33-Rw44+Rw55+Rw66+Rw77+Rw88; 
Rws=Rw15+Rw26-Rw37+Rw48: 
Rwl=Rw13~Rw17—-Rws33~—Rw57+Rw24+Rw28+Rw64+Rw68:.: 
Rwd=Rw14+Rw18—-Rw54-Rw58+Rw23—-Rw27+Rw63+Rw67: 
Rwt=Rw12+Rw1l16-Rw52-Rw56+Rw34+Rw38+Rw744+Rw78: 
Rw=Rwb+Rws+Rw!l-Rwd-Rwt 

Rf=Rfs+Rfp 

Rt=Rw-+Rf 


plot(V.Rt,-".V.Rw,-") 


tO 
U2 
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% ee ee DC EO t** tt tee tte et REESE REESE SHH 


% This file creates the offsets for the pod by obtaining nose 
% tail and parallel midbody information. The nose is formed 
@ using an ellipse and the tail is formed using a parabola. 

% This is then fit to a seventh order polynomial and the results compared. 


seg=input('The min number of segments (must be even) for the pod = '); 
D=input('The max. diameter(ft) of the pod = '); 

L2D=input(‘The desired Legnth to Diameter ratio of the pod = '): 
Lpmb=input(' The amount of Parallel Midbody(ft) in the pod = '): 
Pn=input(' The decimal amount of the pod that is nose section = '): 


Lpod= D*L2D. 
Ln= Pn*(Lpod-Lpmb); 
Lt= Lpod-Ln-Lpmb: 


npn=input(’ The nose shape exponent.(bigger is fuller) =‘); 
npt=input(' The tail shape exponent.(bigger is fuller) = '): 


ppod=zeros( 1,8): 

xn=zeros( |.seg*2):; 
offn=zeros( | .seg*2): 
xpmb=zeros( | .seg/2): 

offpm b=zeros(1.seg/2): 
xt=zeros( 1 .seg*2); 
offt=zeros( 1 .seg*2): 
xpod=zeros( | .(seg*4+seg/2)); 
offpod=zeros(1,(seg*4+seg/2)); 
Xpod=zeros(1,sez*6); 

OFF pod=zeros(1,seg*6): 


xn= (linspace(0,-Ln.seg*2))': 
offn= D/2*(1-((Ln-abs(xn))./Ln).’ npn).“(1/npn): 


xpmb= (linspace(-Ln,-(Ln+Lpmb).seg/2))’; 
offpmb= D/2*xpmb./xpmb.: 


xt = (linspace(-(Ln+Lpmb),-Lpod.seg*2))'; 
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offt= D/2*(1-((abs(xt)-Ln-Lpmb)./Lt).“npt): 


xpod=[xn:xpmb:xt]: 
offpod=[offn:offpmb:offt]; 


ppod=polyfit(xpod.offpod,7)' 


Xpod=(linspace(0,-Lpod.seg*6)): 

OFF pod=(ppod(1).* Xpod.*7)+(ppod(2).* Xpod.*6)+(ppod(3).* Xpod.*5)+... 
(ppod(4).* Xpod.*4)+(ppod(5).* Xpod. 3)+(ppod(6).* Xpod.%2)+... 
(ppod(7).* Xpod.*1)+ppod(8): 


plot(Xpod.OFF pod.’-’.xpod.-OFF pod.’:’) 


24) 


% EEK EK EK EEERES*EECTR LT GE OFS Ee EEE EEE EE EE ERE EEE EEE EE 
% 

% This file creates the offsets for the strut by obtaining nose 

% tail and parallel midbody information. The nose and tail are formed 
% using a parabola. This is then fit to a seventh order polynomial 

% and the results compared. 


segs=input('The min number of segments (must be even) for the strut = '); 
B=input(‘The max. beam/(ft) of the strut = '); 

L2Ds=input(‘The desired Legnth to Diameter ratio of the strut = '): 
Lpmbs=input(' The amount of Parallel Midbody(ft) in the strut = '); 
Pns=input(’ The decimal amount of the strut that is nose section = '): 


Lstrut= B*L2Ds: 
Lns= Pns*(Lstrut-Lpmbs); 
Lts= Lstrut-Lns-Lpmbs; 


nsn=input(' The nose shape exponent,(bigger 1s fuller) = '): 
nst=input(’ The tail shape exponent,(bigger is fuller) = '): 


xns=zeros( | ,segs* 2): 
offns=zeros(1,segs* 2): 
xpmbs=zeros(1.segs/2); 
offpmbs=zeros(1.segs/2): 
xts=zeros(1,segs* 2): 
offts=zeros(1.segs* 2); 
xstrut=zeros(1.(segs*4+segs/2)): 
offstrut=zeros(1.(segs*4+segs/2)); 
Xstrut=zeros(1,segs*6); 

OF Fstrut=zeros( 1 ,segs*6): 


xns= (linspace(0.-Lns,segs*2))’; 
offns= B/2*(1-((Lns-abs(xns))./Lns).“nsn); 


xpmbs= (linspace(-Lns,-(Lns+Lpmbs),segs/2))'; 
offpmbs= B/2*xpmbs./xpmbs: 


xts = (linspace(-(Lns+Lpmbs),-Lstrut,segs*2))’; 


to 
> 
N 


offts= B/2*(1-((abs(xts)-Lns-Lpmbs)./Lts).“nst); 


xstrut=[xns:xpmbs:xts]: 
offstrut=[offns:offpmbs:offts]; 


ps=polyfit(xstrut.offstrut.7 )’ 


Xstrut=(linspace(0.-Lstrut.segs*6)); 

OFF strut=(ps(1).* Xstrut.” 7)+(ps(2).* Xstrut.%6)+(ps(3).* Xstmut.%5)+... 
(ps(4).* Xstrut.“4)+(ps(5).* Xstrut.°3)+(ps(6).* Xstrut.“2)-... 
(ps(7).* Xstrut)—(psi8)): 


plot( Xstrut,OFF strut.'--’ xstrut.offstrut.’-'): 
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oe ee ee ONE ONE m**** ** * * * * * ¥ * & FF EF FF 


thetal=zeros(tseg+1,1!); 
z1=zeros(zseg+1,1); 

x l=zeros(xsegt+1.1); 
plx=zeros(xseg+1,1); 

ql x=zeros(xseg+1,1); 
plz=zeros( zseg+1,1): 
qlz=zeros(zsegt+1.1); 
plt=zeros(tseg+1,1): 
q)t=zeros(tseg+ 1,1); 

res] 1=zeros(tseg+1,1); 
ql At=zeros(tseg+!.1): 

pl At=zeros(tseg+].1); 

pl Az=zeros( zseg+ 1.1): 
ql Az=zeros(zseg+ 1,1); 
z1 A=zeros(zsegz+1,1); 

pl Ax=zeros(xseg+1,1); 
ql Ax=zeros(xseg+1,1).: 

x | A=zeros(xseg+1,1): 
sara]=zeros(xseg+1,1): 
katep 1=zeros(xseg+1].1); 
kateq1=zeros(xsez+1,1); 
greg ]=zeros(xseg+1,1): 
saral A=zeros(xseg+1.1): 
katep] A=zeros(xseg+1.1): 
kateq 1 A=zeros(xseg~+1.1): 
greg | A=zeros(xseg+1,1); 
V=zeros(hs~—1-ls,1); 
Rw11=zerosths+1-ls,1); 


for iv=l:hs: 
V(iv)=Is—(v* lant); 
K=O (is) 2): 
w=-V(1v) (2* pl); 


for it=]:tseg+1-: 
thetal (it)=(pi/2)*((it-1 )/tseg): 


for 1zl= |:zseg+1- 
z1(izl) = -TS*((izl-1)/zseg), 


for 1xl=]:xseg+]; 
x1(x1) =-LS*((ix1-1)/xseg ): 


saral (1x1 )=w*((cs6*x1(1x1)%6)+(cs5*x 1(1x1)°5)+(cs4* x} (1x 1)*%4)+(es3*x 1(1x1)%3)+(es2 
>< (inl 2) (cs l= ix ieeso ie 

katep1(1x1)=cos(k*(sb1 *cos(thetal (it))+os1 *sin(thetal (it)))*(sec(thetal (it)))*2): 

kateq] (1x1 )=sin(k*(sb1 *cos(thetal (it))—os] *sin(thetal (it)))*(sec(theta 1 (it)))* 2): 

oreg 1(1x1)=exp(k* sub 1 *(sec(theta] (it))) *2): 


pl xx] )=saral(ix1)*katep](ix})* greg] (ix} ): 
qi x(ix1)=saral(ix1)*kateq!](ix1)* greg 1(ix1): 


end °othis end closes tne’ **** for 1xl-1_xSeo-—ioopeaoe 


plz(izl )=sum(p1x)/xseg; 
qlz(izl)=sum(q1!x)/xseg: 


end “othis end closes the ~*~" for 1Zl— Zses— Is leopmaauae 


plt(it)=sum(p1z)/zseg: 
ql tit)=sum(qlz)/zseg: 


0% *****e**This starts the next body calculations mside the inctaleop. a am 
for 1z1A= |:zseg—1: 
zlA(izlA) = -TS*((izlA-1)/zseg): 


for 1ixlA=1:xseg+l: 
x]AQx1lA) =-LS*((ix1A-1)/xseg); 


saral A(ix] A)=w*((cs6*x TA(ix1A)*6)+(cs5*x1AGxX1A)°5)-(cs4*#x 1 AGx1lA)y4)+(cs3*x 
1AQx1A)*%3)+(cs2*x1TACx1A)°2)+(ces1 *¥x] ACGix1A))+ces0): 

katep] Aix] A)=cos(k*(sb1*cos(thetal (it))+os1 * sin(thetal (it)))*(sec(thetal (it)))*2); 
kateq] AQx1A)=sin(k*(sb1 *cos(thetal (it))+os1 *sin(thetal (it)))*(sec(theta 1 (it)))*2): 
greg 1 A(ix]A)=exp(k* sub] *(sec(thetal (it)))*2): 


plAx(ixlA)=saral A(x] A)*katep] AGx]A)* greg] A(ix 1A); 
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ql Ax(ix1A)=saral A(ix1A)*kateql] A(ix1A)* greg] A(ix1A); 
end “this end closes the **** for ix] A=1:xseg+1 loop **** 


pl Az(izlA)=sum(p!1Ax)/xseg; 
ql Az(iz1 A)=sum(ql Ax)/xseg; 


end “this end closes the **** for iz1A=1:zseg~1 loop **** 


plAt(it)=sum(p1 Az)/zseg; 
ql At(it)=sum(ql Az)/zseg; 


res1 1 (it)=16* pi*rho*k*2*((plt(it)*pl Atait))+(q1 tiit)*ql] Attit)))*(sec(thetal (it))%3);: 


>) 


end “this end closes the **** for it= l:tseg+1 loop **** 


Rwll(iv)= sum(res1 J )/tseg 


eng oils closes the ~*** foriv = ... loop ***** 
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UF FEE EEREEEEESESHESEECT XN SIX 7 0 ee 


Rw66=zeros(hs-1-ls.1): 
res66=zeros(tseg—1.1): 
thetal=zeros(tseg~1,1): 
zl=zeros(zseg+1.1): 

x 1=zeros(xseg~+1.1): 

pl x=zeros(xseg-1.1): 
qix=zeros(xseg~1.1): 
plz=zeros(zseg~1.1); 

ql z=zeros(zseg~1.1); 
plt=zeros(tseg+1.1): 
qlt=zeros(tseg+1.1): 

ql At=zeros(tseg—1.1): 

pl At=zeros(tseg~1.1): 
pIAz=zeres(zsee— alae 
ql Az=zeros(zseg-1.1): 
zl A=zeros(zseg~ 1.1): 
plAx=zeros(xseg~1.1): 
ql Ax=zeros(xseg—1.1): 

x | A=zeros(xseg-1.1): 
saral=zeros(xseg-~1.1): 
katep 1=zeros(xseg-1.1): 
kateq|1=zeros(xseg—1.1): 
greg l=zeros(xseg-1,1): 
sara l A=zeros(xseg-—1.1): 
katep] A=zeros(xseg~+1.1): 
kateq] A=zeros(xseg~1.1): 
greg | A=zeros(xseg-1,1): 
V=zeros(hs+1-ls.1): 


for iv=1:hs; 
V(iv)=Is+(iv*lani): 
k=g/(V(iv)*2); 
=v a) 
for i1t=1:tseg+]: 
theta] (it)=(pi/2)* ((it-1)/tseg); 


for 1z1= 1:zseg+1: 
z1(iz1) = -TP*((iz1-1)/zseg); 


for ixl=1:xsegtl: 
x1Qx1) =-LP*((ix1-1)/xseg); 


sara] (1x 1 =w*((cp6*x1(1x1)°6)+(cp5*x 1 (1x 1 )*5)+(cp4*x 1(ix 1)*4)+(ep3*x 1x 1)%3)+(ep 
Eoehoix ly 2 )=(epl*xlaxl))rcpd): 

katep | (1x 1 )=cos(k*(sb6*cos(thetal (it))+os6*sin(theta 1 (it)))*(sec(thetal (it)))*2); 

kateq] (1x 1 )=sin(k*(sb6* cos(thetal (it))+os6* sin(theta | (it)))*(sec(thetal (it)))*2); 

greg 1 (1x 1)=exp(k*sub6*(sec( thetal (it)))*2); 


pl xix 1 )=saral (ix1)*katep](ix1)* greg1(ix1): 
ql x(x 1 )=saral(ix1)*kateg!(ix1)* greg 1(ix 1): 


enacts end closes the **** for ixl=1]:xseg~1 loop **** 


plz(izl )=sum(p1x)/xseg: 
qlz(izl )=sum(q|x)/xseg: 


end %this end closes the **** for izl=l:zseg+1 loop **** 


plt(it)=sum(p1 z)/zseg: 
ql t(it)=sum(qlz)/zseg; 


% *******This starts the next body calculations inside the theta loop***** 
for 1z1 A= |:zseg+1: 
ZIA(izlA) = -TP*((iz1A-1)/zseg); 


for 1x l|A=I:xseg~1. 
x]AQx1A) =-LP*((ix1A-1)/xseg); 


saral Ax 1A)=w*((cp6*x1lAQx1A)*6)+(cp5*x 1AGx1A)*5)+(cp4*x1AGx1A)*4)+(cp3* 
mix A) 3)+(cp2*xlAGxlA)*2)+(epl*x] AGixlA))+cp0): 

katep] AGix1A)=cos(k*(sb6* cos(thetal (it))+os6*sin(thetal (it)))*(sec(thetal (it)))*2): 
kateq] A(ix1A)=sin(k*(sb6*cos(thetal (it))+os6* sin(thetal (it)))*(sec(thetal (it)))*2); 
greg 1 A(ixlA)=exp(k*sub6*(sec(thetal (it)))*2); 


plAx(ix1A)=saral A(ix1A)*katep]A(ixlA)* greg] A(ixlA); 
ql Ax(ix1A)=saral AQx1A)*kateq] AGx1A)* greg] A(ixlA); 


end %othis end closes the **** for 1x] A=]:xseg~1 loop **** 
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pl Az(izl1A)=sum(p1Ax)/xseg; 
ql AzizlA)=sum(q] Ax)/xseg; 


end “othis end closes the’***~ Tor 1z/A=) zseeninlooon 
pl Atat)=sum(p1Az)/zseg; 
ql! AtUit)=sum(q1 Az)/zseg; 
res66(it)=16* pi*rho*k*%2*((p1ttit)* pl Atcit))—(q1 tiit)* ql Atat)))*(sec(thetal (it))%3); 


end “othis end closes the *>2* fon it= Miseea ls Gopmaras 


Rw66(iv)= sum(res66)/tseg 


end °othis closes the **** {ori ==, loop. 7am 
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Dee NIE SX W474 tt tt tt tt tt tt te ee eee 


Rw16=zeros(hs+1-ls,1); 
res 16=zeros(tseg+1,1); 
thetal=zeros(tseg+1,1); 
zl=zeros(zseg+1,1); 

x ]=zeros(xseg+1,1); 
p|x=zeros(xseg+1,1); 

ql x=zeros(xseg+1,1); 
plz=zeros(zseg+1,1); 
qiz=zeros(zseg+1,1); 
pilt=zeros(tseg+1.1\); 
qlt=zeros(tseg+ 1,1); 

ql At=zeros(tseg+1,1); 
plAt=zeros(tseg+1,1); 

pl Az=zeros(zseg+1,1); 
ql Az=zeros(zseg+ | ,1): 
zl A=zeros(zseg+1.1); 
plAx=zeros(xseg+ 1.1): 
ql Ax=zeros(xseg+1,1); 

x | A=zeros(xseg+1,1); 
sara | =zeros(xseg+1,1); 
katep 1=zeros(xseg+1.1); 
kateql=zeros(xseg+ 1,1); 
greg |=zeros(xseg+1,1); 
saral A=zeros(xseg+1,1): 
katep] A=zeros(xseg+1,1); 
kateq] A=zeros(xseg+1,1); 
greg|1 A=zeros(xseg+1,1): 
V=zeros(hs+1-ls.1): 


for 1iv=1:hs; 
V(iv)=Is~(iv*lant); 
k=g/(V(iv)*2); 
=-V(iv)/(2* pi): 
for 1t=1:tseg+1; 
thetal (it)=(pi/2)*((it-1)/tseg); 


for izl= |:zseg+1: 
z1(iz1) = -TP*((1iz1-1)/zseg): 
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10 1x) —lexseoael 
x1l(ixl) =-LP*((ix1-1)/xseg); 


saral (ix 1 )=w*((cs6*x 1 (1x 1)*6)+(cs5*x1(1x1)*5)+(cs4*x 1(1x1)*4)+(cs3*x 1 (1x 1)%3)+(es2 
*x 10x] )*2)+(csl*x1Gxt))rcs0): 

katep1(1x1 )=cos(k*(sb1 *cos(thetal (it))+os1 * sin(thetal (it)))*(sec(thetal (it)))*2); 

kateq| (1x 1)=sin(k*(sb1 *cos(thetal (1t))+os1 * sin(thetal (it)))*(sec(thetal (it)))*2); 

greg] (ixl )=exp(k* sub] *(sec(thetal (1t)))*2); 


plx(ix1)=saral(ix1)*katep](ix1)* greg] (ix); 
qi x(ix})=saral (ix1)*kateq](ix1)* greg1(ix1); 


end “this end closes the **** for ixl=1l:xseg+1 loop **** 


plz(izl)=sum(p1x)/xseg; 
qi z(iz1)=sum(q!x)/xseg; 


end “othis end closes the **** for izl1=1:zseg+1 loop **** 


plt(it)=sum(p1z)/zseg; 
qlit(it)=sum(q]z)/zseg: 


% *******T his starts the next body calculations inside the theta loop***** 
for 1z1A= 1|:zseg+1; 
z1A(QzlA) = -TP*((izlA-1)/zseg); 


for 1x1 A=1]:xsegt+1]; 
x1A(ixlA) =-LP*((ix1A-1)/xseg); 


saral A(ix]A)=w*((cp6*x1A(ix1A)*%6)+(cp5*x1A(ixlA)*5)+(cp4*x1lA(ix 1A)*4)+(cp3* 
x] A(ixlA)%3)+(cp2*xlA(UixlA)*2)+(cp1*¥x 1] A(QixlA))+cp0); 
katepl1A(ix1A)=cos(k*(sb6*cos(thetal (it))+0s6* sin(thetal (it)))*(sec(thetal (1t)))*2); 
kateq] A(ix] A)=sin(k*(sb6*cos(theta | (it))+0s6* sin(thetal (it)))*(sec(theta 1 (1t)))*2); 
greg] A(ixl A)=exp(k*sub6*(sec(thetal (it)))*2); 


plAx(ixl]A)=saral AGx1A)*katep]A(ixlA)* greg] Aix] A); 
qlAx(ixlA)=saral AGix1A)*kateq] Aix] A)* greg] A(ixlA); 


end “%this end closes the **** for ix] A=1:xseg+! loop **** 
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plAz(iz1A)=sum(p] Ax)/xseg; 
ql Az(izl A)=sum(q!Ax)/xseg; 


end “this end closes the **** for 1z1A=1:zseg+1 loop **** 


pl At(it)=sum(p1Az)/zseg: 
qi At(it)=sum(q1Az)/zseg: 


res16(it)=32* pi*rho*k*2*((pltiit)* pl At(it))+(q]1 t(it)* ql At(it)))* (sec(thetal (it))%3); 


end “this end closes the **** for it= 1:tseg+1 loop **** 


Rw16(iv)= sum(res!6)/tseg 


em@e, -o this closes the **** for iv =... loop ***** 
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